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ABSTRACT 

We present an updated description of the Planck Low Rrequency Instrument (LRl) data processing pipeline, associated with the 2015 data release. 
We point out the places where our results and methods have remained unchanged since the 2013 paper and we highlight the changes made for 
the 2015 release, describing the products (especially timelines) and the ways in which they were obtained. We demonstrate that the pipeline is 
self-consistent (principally based on simulations) and report all null tests. Ror the first time, we present LRI maps in Stokes Q and U polarization. 
We refer to other related papers where more detailed descriptions of the LRI data processing pipeline may be found if needed. 

Key words. Space vehicles: instruments - Methods:data analysis - cosmic microwave background 


1. Introduction 

This paper, one of a set associated with the 2015 release of 
data from the Planck^ mission (Planck Collaboration I 2016), 
describes the Low Frequency Instrument (LFI) data processing 
that supports the second Planck cosmological release. Following 
the nominal mission of 15.5 months, the LFI in-flight operation 
was extended to fully exploit the lifetime of the Planck 20 K to 


* Corresponding author: A. Zacchei, zacchei@oats. inaf. it 
^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from Prance and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


4 K cryogenic system, leading to a total of 48 months of ob¬ 
servation (or eight full-sky surveys) with essentially unchanged 
instrument performance. This paper is an updated description 
of the LFI data processing (Planck Collaboration II 2014) that 
was part of the second wave of astrophysical results published 
in early 2014 (Planck Collaboration VIII-XXVI 2013), now in¬ 
corporating the analysis of the full mission data, both in temper¬ 
ature and in polarization. This work describes the overall data 
flow of the pipeline implemented at the LFI data processing cen¬ 
tre (DPC), including scientiflc telemetry from the instrument, 
housekeeping data, and frequency maps, as well as the tests ap¬ 
plied to validate the data products. Detailed descriptions of criti¬ 
cal aspects of the data analysis and products, including improve¬ 
ments in some of the algorithms used in the pipeline, are given 
in four companion papers. These discuss, respectively: system¬ 
atic effects and the overall error budget (Planck Collaboration III 
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2016); the determination of the LFI main beams and window 
functions from in-flight planet-crossing measurements and opti¬ 
cal modelling (Planck Collaboration IV 2016); photometric cal¬ 
ibration, including methods adopted and related uncertainties 
(Planck Collaboration V 2016); and mapmaking, including the 
process used to obtain the low-resolution maps and their asso¬ 
ciated full noise covariance matrices (Planck Collaboration VI 
2016). The main results and reference tables on all these topics 
are summarized in this paper. 

The structure of this paper is organized as follows. We sum¬ 
marize the overall data processing pipeline in Sect. 3. Processing 
of the time ordered information (TOI) is described in Sect. 4, 
with an emphasis on changes since Planck Collaboration II 
(2014). Section 6 describes important changes to our calcula¬ 
tions of LFI beams, which in turn has an effect on calibra¬ 
tion, described in Sect. 7. LFI noise properties are described in 
Sect. 8. Section 9 and Sect. 10 present Planck maps at 30, 44, 
and 70 GHz, both in temperature and in Q and U polarization, 
including the low-multipole maps needed to construct the Planck 
likelihood (Planck Collaboration XI 2016). Section 11 presents 
the major new results for this release, LFI polarization maps, 
and an analysis of systematic effects peculiar to polarization. 
Validation of the LFI products, especially by means of null tests, 
is discussed in Sect. 12, and the special issue of data selection for 
low-^ analysis is considered in Sect. 13. Section 15 summarizes 
the LFI data products (for further details, see the Explanatory 
Supplement^ that accompanies the release of products and pro¬ 
vides its detailed description). We conclude briefly in Sect. 16. 

2. In-flight behaviour and operations 

The Planck LFI instrument is described in Bersanelli et al. 
(2010) and Mennella et al. (2010). It comprises 11 radiome¬ 
ter chain assemblies (RCAs), two at 30 GHz, three at 44 GHz, 
and six at 70 GHz, each composed of two independent pseudo¬ 
correlation radiometers sensitive to orthogonal linear polariza¬ 
tion modes. Each radiometer has two independent square-law 
diodes for detection, integration, and conversion from radio fre¬ 
quency signals to DC voltages. The focal plane is cryogenically 
cooled to 20 K, while the pseudo-correlation design uses inter¬ 
nal, blackbody reference loads cooled to 4.5 K. The radiome¬ 
ter timelines are produced by taking differences between the 
signals from the sky, and from the reference loads, Vr^f. 
Radiometer balance is optimized by introducing a gain modula¬ 
tion factor, typically stable within 0.02 % throughout the entire 
mission, which greatly reduces 1 // noise and improves immu¬ 
nity to a wide class of systematic effects (Mennella et al. 2011). 
During the full operation period (ignoring a brief, less stable 
thermal period due to the sorption cooler switchover), the be¬ 
haviour of all 22 LFI radiometers was stable, with 1 // knee fre¬ 
quencies unchanging within 9 % and white noise levels within 
0.5 %. These results are in line with those found for the 15.5 
month nominal mission period (Planck Collaboration II 2014). 

2.1. Operations 

The data set released together with this paper was acquired from 
12 August 2009 to 3 August 2013, roughly four years of obser¬ 
vations. The first two years of data (from Survey 1 to Survey 4) 
were acquired scanning the sky with a phase angle of 340°, 
whereas for the last two years (from Survey 5 to Survey 8) 
the phase angle was shifted to 250° (see Planck Collaboration I 

^ http://pla.esac.esa.int/pla/index.html 


2016, for details). This shift has allowed for more thorough in¬ 
vestigation of systematic effects, including better characteriza¬ 
tion of the beam and the related Galactic straylight (see Sect. 7.4) 
using null tests based on survey differences. During the last three 
Jupiter crossing, the scanning strategy was optimized to obtain 
a better beam determination (see Sect. 6). The period from 03 
August 2013 to 03 October 2013 was used to perform deep scan¬ 
ning of the Crab Nebula and of the regions near the minima of 
the cosmic microwave background (CMB) dipole, with the aim 
of determining the dipole direction with an alternative approach. 
These data are not included in this release, since they require 
specialized analysis, which is not yet complete. 

2.2. Instrument performance update 

In Table 1 we present a top-level summary of the instrument 
performance parameters measured in flight during the four 
years of operation of LET Optical properties have been recon¬ 
structed from Jupiter transits (Planck Collaboration IV 2016) 
and are in agreement with estimations made for the 2013 re¬ 
lease (Planck Collaboration IV 2014). White noise sensitivity 
and parameters describing the 1// noise component are in line 
with the 2013 values (Planck Collaboration II 2014), demon¬ 
strating that cryogenic operation of the low-noise amplifiers and 
phase switches do not result in any significant aging effects over 
a period of four years. Overall calibration uncertainty, deter¬ 
mined as the sum of absolute and relative calibration, is 0.35 %, 
0.26 %, and 0.20 % at 30, 44, and 70 GHz respectively, improv¬ 
ing by more than a factor of 2 over the LFI 2013 calibration 
(Planck Collaboration V 2014). The residual systematic uncer¬ 
tainty was computed for both temperature and polarization; it 
varies between 1 and 3yuKcMB (Planck Collaboration III 2016) 
in temperature and polarization. It should be noted that the un¬ 
certainty arising from systematic effects is lower than in the pre¬ 
vious release (Planck Collaboration IV 2014); this is principally 
due to the straylight removal and the new iterative calibration 
algorithm now used. 

3. Data processing overview 

As in Planck Collaboration II (2014), the processing of LFI data 
is divided into three levels, shown schematically in Fig. 1. The 
main changes compared to the earlier release are related to the 
way in which we take into account the beam information in the 
pipeline processing, as well as an entire overhaul of the iterative 
algorithm used to calibrate the raw data. According to the LFI 
scheme, processing starts at Level 1, which retrieves all the nec¬ 
essary information from data packets and auxiliary data received 
from the Mission Operation Centre, and transforms the scien¬ 
tific packets and housekeeping data into a form manageable by 
Level 2. Level 2 uses scientific and housekeeping information 
to: 

- build the LFI reduced instrument model (RIMO), which con¬ 
tains the main characteristics of the instrument; 

- remove analogue-to-digital converter (ADC) non-linearities 
and 1 Hz spikes diode by diode (see Sects. 4.2 and 4.3); 

- compute and apply the gain modulation factor to minimize 
1 If noise (see Sect. 4.4); 

- combine signals from the two diodes of each radiometer (see 
Sect. 4.5); 

- compute the appropriate detector pointing for each sample, 
based on auxiliary data and beam information corrected by 
a model (PTCOR) built using solar distance and radiometer 
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Table 1. LFI performance parameters. 


Parameter 

30 GHz 

44 GHz 

70 GHz 

Centre frequency [GHz]. 

28.4 

44.1 

70.4 

Scanning beam FWHM^ [arcmin]. 

33.10 

27.94 

13.08 

Scanning beam ellipticity^ . 

1.37 

1.25 

1.27 

Effective beam FWHM^ [arcmin] . 

32.29 

27.00 

13.21 

White noise level in timelines^ [//Kcmb s^^^]. 

148.1 

174.2 

152.0 

fkn^ [mHz]. 

114 

52 

19 

1 If slope^. 

-0.92 

-0.88 

-1.20 

Overall calibration uncertainty^ [%]. 

0.35 

0.26 

0.20 

Systematic effects uncertainty in Stokes P [^uKcmb] • 

0.88 

1.97 

1.87 

Systematic effects uncertainty in Stokes g® [)uKcmb] 

1.11 

1.14 

2.25 

Systematic effects uncertainty in Stokes [)uKcmb] 

0.95 

1.20 

2.22 


^ Determined by fitting Jupiter observations directly in the timelines. 

^ Calculated from the main beam solid angle of the effective beam (Sect. 6.2). These values are used in the source extraction 
pipeline (Planck Collaboration XXVI 2016). 

^ Typical values derived from fitting noise spectra (Sect. 8.1). 

^ Sum of the error determined from the absolute and relative calibration, see Planck Collaboration IV (2016). 

^ Peak-to-peak difference between 99 % and 1 % quantiles in the pixel value distributions from simulated maps (see Planck Collaboration III 
2016). 


electronics box assembly (REBA) temperature information 
(see Sect. 5); 

“ calibrate the scientific timelines to physical units (Kcmb), fit¬ 
ting the total CMB dipole convolved with the Ttt beam rep¬ 
resentation (see Sect. 7), without taking into account the sig¬ 
nature due to Galactic stray light (see Sect. 7.4); 

- remove the solar and orbital dipole convolved with the Ttt 
beam representation and the Galactic emission convolved 
with the beam sidelobes (see Sect. 7.4) from the scientific 
calibrated timeline; 

- combine the calibrated time-ordered information (TOI) into 
aggregate products, such as maps at each frequency (see 
Sect. 9). 

Level 3 collects Level 2 outputs from both HFI 
(Planck Collaboration VIII 2016) and LFI and derives various 
products, such as component-separated maps of astrophysical 
foregrounds, catalogues of different classes of sources, and the 
likelihood of cosmological and astrophysical models given in 
the maps. 

4. Time-ordered information (TOI) processing 

The Level 1 pipeline, which has the responsibility to receive 
telemetry data and sort them into a form manageable by the 
Level 2 pipeline, has not changed with respect to the 2013 re¬ 
lease; we therefore refer to Planck Collaboration II (2014) for 
its description. In this section, we move directly to a discussion 
of the Level 2 pipeline. 

4.1. Input flags 

The flagging procedure used was exactly the same as de¬ 
scribed in Planck Collaboration II (2014). In Table 2 we give 
the percentage of usable and unused data for the full mis¬ 
sion. It should be noted that compared with the same table 
in Planck Collaboration II (2014) the amount of missing data 
(where by “missing” we mean packets that were not been re¬ 
ceived on the ground) is larger due to two technical problems 


Table 2. Percentage of LFI observation time lost due to missing 
or unusable data, and to manoeuvres^. 


Category 

30 GHz 

44 GHz 

70 GHz 

Missing [%]. 

0.153 

0.154 

0.153 

Anomalies [%]. 

0.375 

0.448 

0.631 

Manoeuvres [%] . . . . 

8.032 

8.032 

8.032 

Usable [%] . 

91.440 

91.366 

91.184 


The remaining percentage (listed in the last row) is used in scientific 
analysis. 


that were experienced with the spacecraft, resulting in data not 
being downloaded for 2 days of observation. On the other hand 
the anomalies were lessened due to better control of the instru¬ 
ment’s temperature stability. The percentages of time spent on 
spacecraft manoeuvres are the same for the three frequencies, 
and as a consequence the fraction of data used in the science 
analysis was similar (at more than 90 %) at each frequency. 

4.2. ADC non-linearity correction 

The ADCs convert the analogue detector voltages to numbers, 
their linearity is as important as that of the receivers and de¬ 
tectors, with any departure appearing as a distortion in the 
system power response curve. While the algorithm for deter¬ 
mining the ADC corrections remains the same as described in 
Planck Collaboration IV (2014), some changes were made in its 
implementation and execution. First, the full mission data are 
now used, so that when detector voltages are revisited there will 
be an improvement in signal to noise (although in the partic¬ 
ular case of radiometer 21M, some of the voltages were too 
poorly sampled to generate an adequate solution). Second, in¬ 
stead of determining the white noise amplitude via a Fourier 
transform, we now use the difference between the sum of the 
variances and twice the covariance of adjacent paired points in 
the timestream, such that white noise variance = Var[Xo] + 
Var[Xe] - 2Cov[Xo,Xe], where Xq and Xq are data points with 
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Fig. 1. Schematic representation of the Level 2 and pointing pipelines of the LFI DPC; elements in red identify those modified or 
augmented with respect to Planck Collaboration II (2014) 


odd and even indices, respectively. This not only increased the 
speed of calculating the noise amplitude, but avoided the itera¬ 
tion steps, since these can be done analytically from the initial 
variance-covariance estimates. Finally, data acquisition electron¬ 
ics (DAE) offset changes made on operational day (OD) 953 to 
avoid saturation also shifted the apparent ADC voltage relative 
to the true detector voltage. A separate ADC correction had to 
be generated and applied to radiometers 22M and 23 S using only 
the post OD 953 data. 


The ability to recover the correct ADC solution and the level 
of the residuals was assessed by simulating time-ordered data 
with the same noise statistics, voltage drift, gain fluctuations, 
and sky signal, with a known ADC error. As the correction in the 
DPC pipeline is a lookup table of input to output detector volt¬ 
ages to which a spline is fitted to interpolate the TOI voltages, we 
introduced the ADC error as the spline curve with the input and 
output voltages swapped and thereby generate the inverse of the 
measured ADC effect. Comparing the spline curves used to the 
ones recovered proved to be at the level of a few percent, lead¬ 
ing to rms errors on the residual simulated frequency maps of 
0.1 yuKcMB at 30 and 44 GHz and 0.4 yuKcMB at 70 GHz, both 
temperature and polarization. These simulations and results are 
summarized in more detail in Planck Collaboration III (2016). 


4.3. Corrections for electronic spikes 

Electronic spikes are caused by the interaction between the elec¬ 
tronics clock and the scientific data lines. They occur in the 
data acquisition electronics (DAE) after the detector diodes and 
before the analogue-to-digital converters (ADC, Meinhold et al. 
2009; Mennella et al. 2010, 2011). The signal is detected in all 
the LEI radiometers time-domain outputs as a 1 s square wave 
with a rising edge near 0.5 s and a falling edge near 0.75 s, syn¬ 
chronous with the on-board time signal. In the frequency domain 
it appears as a spike signal at multiples of 1 Hz. The 44 GHz 
channels are the only one that are significantly affected by this 
effect. Consequently the spike signal is removed from the data 
only in this channel. The procedure consists of the subtraction 
of a fitted square wave template from the time-domain data as 
described in Planck Collaboration III (2016). We are evaluating 
the possibility of further reducing the residual effect of the spikes 
signal at the map level, as described in Planck Collaboration III 
(2016), for the next Planck data release by adopting one or more 
of the following approaches: 

- increasing the resolution of square wave template, at the mo¬ 
ment at 80 Hz; 

- using time varying template instead of the fixed one over the 
whole mission; 


4 




































































































Planck Collaboration: LFI data processing 


Table 3. Weights used in combining diodes^. 


Diode 


Radiometer 

M-00 

M-01 

S-10 

S-11 

70 GHz 

LFI 18 

0.567 

0.433 

0.387 

0.613 

LFI 19 

0.502 

0.498 

0.551 

0.449 

LFI 20 

0.523 

0.477 

0.477 

0.523 

LFI 21 

0.500 

0.500 

0.564 

0.436 

LFI 22 

0.536 

0.464 

0.554 

0.446 

LFI 23 

0.508 

0.492 

0.362 

0.638 

44 GHz 

LFI 24 

0.602 

0.398 

0.456 

0.544 

LFI 25 

0.482 

0.518 

0.370 

0.630 

LFI 26 

0.593 

0.407 

0.424 

0.576 

30 GHz 

LFI 27 

0.520 

0.480 

0.485 

0.515 

LFI 28 

0.553 

0.447 

0.468 

0.532 


A perfect instrument would have weights of 0.500 for both diodes. 


- removing spikes signal from the 30 GHz and 70 GHz chan¬ 
nels. 

4.4. Demodulation: gain modulation factor estimation and 
application 

Each Planck LFI diode switches at 4096 Hz (Mennella et al. 
2010) between the sky and the 4 K reference load. The data ac¬ 
quired in this way are dominated by 1/f noise that is highly 
correlated between the two streams (Bersanelli et al. 2010); dif¬ 
ferencing those streams results in a strong reduction of 1// 
noise. The procedure applied differs from that discussed in 
Planck Collaboration II (2014) in only one way: the Galaxy and 
point sources are masked from the time-ordered data used in the 
computation of the gain modulation factor R (GMF in Fig. 1). 
The overall variation of R over the whole mission is less than 
0.02 % for every LFI channel. A full description of the theory of 
this correction can be found in Mennella et al. (2011). 

4.5. Combining diodes 

Two detector diodes provide the output for each Planck LFI 
receiver channel. To minimize the impact of imperfect isola¬ 
tion between the two diodes, we perform a weighted average 
of the time-ordered data from the two diodes of each receiver 
just before the differencing. The procedure applied is the same 
described in Planck Collaboration II (2014) and for the sake of 
completeness, we report in Table 3 the values of weights used; 
the receiver channels are indicated either with M (main) or S 
(side). The weights are kept fixed during the whole mission. 


5. Pointing 

The long time scale pointing correction, PTCOR, has been mod¬ 
ified, and is now based on the solar distance and radiometer 
electronics box assembly (REBA) thermometry. Unlike in 2013, 
the reconstructed satellite attitude is now uniform across both of 
the Planck instruments and is discussed in detail in the mission 
overview paper, Planck Collaboration I (2016). 


6. Main beams and the geometrical calibration of 
the focal plane 

The in-flight assessment of the LFI main beams relies on the 
measurements performed during seven Jupiter crossings; the 
first four transits (“Jl” to “J4”) occurred in nominal scan mode 
(spin shift 2 arcmin, 1 deg per day), and the last three scans 
(“J5” to “J7”) in a deeper coverage mode (spin shift 0.5 arcmin, 
15 arcmin per day). The period of time corresponding to each 
Jupiter observation is reported in Table 4. By stacking data from 
the seven scans, we measure the main beam profiles down to 
-25 dB at 30 and 44 GHz, and down to -30 dB at 70 GHz. If we 
fit the main beam shapes with elliptical Gaussian profiles, the un¬ 
certainties of the measured scanning beams can be expressed in 
terms of statistical errors on these Gaussian parameters. With re¬ 
spect to the 2013 release, the improvement in the signal-to-noise 
ratio due to the number of samples and to better sky coverage is 
about a factor of 2. The beam full width half maximum is de¬ 
termined with a typical uncertainty of 0.2 % at 30 and 44 GHz, 
and 0.1 % at 70 GHz, approximately a factor of 2 better than 
the value achieved in 2013. The fitting procedure also returns 
the main beam pointing directions in the Planck field of view 
(i.e. the focal plane geometry), centred along the nominal line of 
sight as defined in Tauber et al. (2010). 

We determined the focal plane geometry of LFI indepen¬ 
dently for each Jupiter crossing (Planck Collaboration IV 2016), 
using the same procedure as adopted in the 2013 release. The so¬ 
lutions for the seven crossings agree within 4 arcsec at 70 GHz, 
and 7 arcsec at 30 and 44 GHz. The uncertainty in the deter¬ 
mination of the main beam pointing directions evaluated from 
the single scans is about 4 arcsec for the nominal scans, and 
2.5 arcsec for the deep scans at 70 GHz (27 arcsec for the nom¬ 
inal scan and 19 arcsec for the deep scan, at 30 and 44 GHz). 
Stacking the seven Jupiter transits, the uncertainty in the re¬ 
constructed main beam pointing directions becomes 0.6 arcsec 
at 70 GHz and 2 arcsec at 30 and 44 GHz. With respect to the 
2013 release, we have found a difference in the main beam 
pointing directions of about 5 arcsec in the cross-scan direc¬ 
tion and 0.6 arcsec in the in-scan direction. The beam centres 
and polarization orientation are defined by four parameters, 6uv 
and (puv^ which define the beam pointing reconstructed using the 
stacked Jupiter transits; and ij/nv and defining the polariza¬ 
tion orientation of the beam (see Planck Collaboration IV 2016 
and Planck Collaboration ES 2013 for the definitions of these 
angles); their values for all the LFI radiometers are reported in 
Table 5. Only Ouv and (puv, which are the beam pointing in spher¬ 
ical coordinates referred to the line of sight, can be determined 
with Jupiter observations. The polarization orientation of the 
beams, defined by ij/uv + <Apoh is estimated based on the geometry 
of the waveguide components in the LFI focal plane (which for 
coherent detectors defines the polarization planes to high preci¬ 
sion), reprojected in the sky through our GRASP model. As dis¬ 
cussed in Planck Collaboration III (2016), direct measurements 
of bright polarized sources (such as the Crab Nebula) provide 
only loose constraints, and our final uncertainties on the polar¬ 
ization angles have been evaluated through simulations. 

Details of the LFI main beam reconstruction and focal plane 
geometry evaluation are reported in Planck Collaboration IV 
(2016). 

6.1. Scanning beams 

The “scanning beams”, see Table 6 for main beam descriptive 
parameters, used in the LFI pipeline (affecting calibration, ef- 
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Table 4. Approximate dates of the Jupiter observations. 


Jupiter transit 

Date 

Scan 1 (Jl) 

21 Oct - 5 Nov, 2009 

Scan 2 (J2) 

27 Jun - 12 Jul, 2010 

Scan 3 (J3) 

3-18 Dec, 2010 

Scan 4 (J4) 

30 Jul - 8 Aug, 2011 

Scan 5 (J6) 

8-30 Jan, 2012 

Scan 6 (J6) 

1 - 14 Sept, 2012 

Scan 7 (J7) 

7-28 Feb, 2013 


Table 5. Focal plane geometry. 


Radiometer 

0u/ 


tl'uv'" 



[deg] 

[deg] 

[deg] 

[deg] 

70 GHz 

LFI 18M 

3.334 

-131.828 

22.15 

90.2 

LFI 18S 

3.334 

-131.820 

22.15 

0.0 

LFI 19M 

3.209 

-150.482 

22.40 

90.0 

LFI 19S 

3.209 

-150.488 

22.40 

0.0 

LFI 20M 

3.184 

-168.182 

22.38 

89.9 

LFI 20S 

3.185 

-168.194 

22.38 

0.0 

LFI 21M 

3.186 

169.281 

-22.38 

90.1 

LFI 21S 

3.185 

169.271 

-22.38 

0.0 

LFI 22M 

3.174 

151.360 

-22.34 

90.1 

LFI 22S 

3.174 

151.371 

-22.34 

0.1 

LFI 23M 

3.281 

132.259 

-22.08 

89.7 

LFI 23S 

3.281 

132.280 

-22.08 

0.0 

44 GHz 

LFI 24M 

4.073 

-179.540 

0.01 

90.0 

LFI 24S 

4.071 

-179.505 

0.01 

0.0 

LFI 25M 

4.984 

61.093 

-113.23 

89.5 

LFI 25S 

4.983 

61.125 

-113.23 

0.0 

LFI 26M 

5.036 

-61.670 

113.23 

90.5 

LFI 26S 

5.036 

-61.675 

113.23 

0.0 

30 GHz 

LFI 27M 

4.346 

153.987 

-22.46 

89.7 

LFI 27S 

4.346 

153.985 

-22.46 

0.0 

LFI 28M 

4.376 

-153.424 

22.45 

90.3 

LFI 28S 

4.375 

-153.418 

22.45 

0.0 


^ Beam pointing reconstructed using the stacked Jupiter transits. 
^ Polarization orientation of the beam derived from simulations. 


fective beams, and beam window functions) are very similar 
to those presented in Planck Collaboration IV (2014): they are 
GRASP beams properly smeared to take into account the satellite 
motion. They come from a tuned optical model and represent 
the most realistic fit to the available measurements of the LFI 
main beams. These beams have now been validated using seven 
Jupiter transits. The Jupiter scans allow us to measure the total 
field, that is the co- and cross-polar components combined in 
quadrature. The adopted beam model has the added advantage 
that it allows the co- and cross-polar pattern to be defined sepa¬ 
rately; it also permits us to properly consider the beam cross¬ 
polarization in every step of the LFI pipeline. The scanning 
beams reconstructed from Jupiter transits are shown in Fig. 2 
Unlike in Planck Collaboration IV (2014), where the main 
beams were full-power main beams and the resulting beam win¬ 
dow functions were normalized to unity (because the calibra¬ 
tion was performed assuming a pencil beam), a different beam 
normalization is introduced here to properly take into account 


the power entering the main beam (typically about 99% of 
the total power). Indeed, as described in Planck Collaboration V 
(2016), the current LFI calibration takes into account the full 471 
beam (i.e. the main beam, as well as near and far sidelobes). 
Consequently, in the calculation of the window function, the 
beams are not normalized to unity; instead, their normalization 
takes into account the real efficiency calculated by considering 
the variation across the band of the optical response (coupling 
between feedhorn pattern and telescope) and the radiometric re¬ 
sponse (band shape). This affects flux densities derived from the 
maps (see Sect. 7.2). 

In addition, “hybrid beams” have been created using planet 
measurements above 20 dB from the main beam power peak and 
GRASP beams below this threshold. The hybrid beams have been 
normalized to match the GRASP beams (i.e. the main beam effi¬ 
ciency is set to be the same). Hybrid beams have been used to 
perform a further check on the consistency between the GRASP 
model and the planet data, in terms of window functions. Further 
details are reported in Planck Collaboration IV (2016). 



Fig. 2. Scanning beams reconstructed from Jupiter observations. 
The beams are plotted in logarithmic contours of -3, -10, -20, 
and -30 dB from the peak, for the 70 GHz channel (horns 18- 
23), and -3, -10, -20, and -25 dB from the peak, for the 30 
and 44 GHz channel (horns 27 and 28, and 24-26, respectively). 
The main and side arms are indicated with black and blue lines, 
respectively. 


6.2. Effective beams 

The GRASP combined co- and cross-polar main beams are used 
to calculate the “effective beams”, which take into account 
the specific scanning strategy and pointing information in or¬ 
der to include any smearing and orientation effects on the 
beams themselves. We compute the effective beam at each LFI 
frequency, using the scanning beam and scan history in real 
space using the FEBeCoP (Mitra et al. 2011) method. Effective 
beams are used to calculate the effective beam window func¬ 
tion, as reported in (Planck Collaboration IV 2016) and in the 
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Table 6. Main beam descriptive parameters of the scanning 
beams, with +1 cr uncertainties, represents the beam orien¬ 
tation as defined in Planck Collaboration IV (2014) 


Beam 

FWHM 

[arcmin] 

Ellipticity 

0ell 

[degrees] 

70 GHz 

LFI 18M 

13.40 + 0.02 

1.235 + 0.004 

85.74 + 0.41 

LFI 18S 

13.46 + 0.02 

1.278 + 0.004 

86.41 + 0.33 

LFI 19M 

13.14 + 0.02 

1.249 + 0.003 

78.82 + 0.35 

LFI 19S 

13.09 + 0.02 

1.281 + 0.002 

79.15 + 0.30 

LFI 20M 

12.83 + 0.02 

1.270 + 0.003 

71.59 + 0.32 

LFI 20S 

12.83 + 0.02 

1.289 + 0.004 

72.69 + 0.31 

LFI 21M 

12.75 + 0.02 

1.280 + 0.003 

107.99 + 0.27 

LFI21S 

12.86 + 0.02 

1.294 + 0.003 

106.96 + 0.29 

LFI 22M 

12.92 + 0.02 

1.264 + 0.003 

101.87 + 0.30 

LFI 22S 

12.99 + 0.02 

1.279 + 0.003 

101.61 + 0.30 

LFI 23M 

13.32 + 0.02 

1.235 + 0.004 

93.53 + 0.40 

LFI 23S 

13.33 + 0.02 

1.279 + 0.004 

93.49 + 0.36 

44 GHz 

LFI 24M 

23.18 + 0.05 

1.388 + 0.005 

89.82 + 0.33 

LFI 24S 

23.03 + 0.04 

1.344 + 0.003 

89.97 + 0.34 

LFI 25M 

30.02 + 0.07 

1.191 + 0.005 

115.95 + 0.75 

LFI 25S 

30.79 + 0.07 

1.188 + 0.005 

117.70 + 0.74 

LFI 26M 

30.13 + 0.08 

1.191 + 0.006 

61.89 + 0.84 

LFI 26S 

30.52 + 0.08 

1.189 + 0.006 

61.53 + 0.77 

30 GHz 

LFI 27M 

32.96 + 0.06 

1.364 + 0.005 

101.20 + 0.34 

LFI 27S 

33.16 + 0.07 

1.379 + 0.005 

101.29 + 0.34 

LFI 28M 

33.17 + 0.07 

1.366 + 0.006 

78.17 + 0.36 

LFI 28S 

33.12 + 0.07 

1.367 + 0.005 

78.47 + 0.33 


7. Photometric calibration 

With the term “photometric calibration,” we indicate the process 
that converts the raw voltages V measured by the LFI radiome¬ 
ters into a thermodynamic temperature. The response of an LFI 
radiometer to a change in the temperature coming from the sky 
can be modelled by the following equation: 

y(0 = Gx[B*(£) + 7’cMB+7’sky) + 7’o], (1) 

where B is the beam response, the temperature T = D Tcmb + 
Tsky is decomposed into the sum of three terms (the dipole in¬ 
duced by the motion of the Solar system plus the Planck space¬ 
craft, the CMB, and any other foregrounds), and Tq is a constant 
offset, which includes both instrumental offsets and the CMB 
monopole. The quantity G is the unknown term in the calibration 
problem, and its inverse K = G“\ the “calibration constant,” is 
used to convert the timestream of voltages V(t) into tempera¬ 
tures. 

Planck's calibration source has always been the dipole 
term, D. However, since the previous Planck data release 
(Planck Collaboration V 2014) we have implemented a number 
of important changes in the pipeline used to calibrate the volt¬ 
ages measured by the LFI radiometers. In this section we provide 
an overview of the most important result; we refer the reader in¬ 
terested in further details to Planck Collaboration V (2016). 

We use as a calibration source the signal B ^ D'm Eq. (1), 
which is induced by the combined motion of the spacecraft and 
the Solar System with respect to the CMB rest frame. We have 
characterized the dipole by means of Planck data and have es¬ 
timated the amplitude to be (3364.5 + 2.0)yuKcMB in the direc¬ 
tion I = 264?00 + 0?03, b = 48?24 + 0?20 in Galactic coordi¬ 
nates (Planck Collaboration I 2016). This represents an approxi¬ 
mately 0.3 % increase in the amplitude with respect to the dipole 
used in the 2013 data release, which was based on the results of 
Hinshaw et al. (2009). 


source detection pipeline used to generate the PCCS catalogue 
(Planck Collaboration XXVI 2016). Table 7 lists the mean and 
rms variation across the sky of the main parameters computed 
with FEBeCoP. Note that the FWHM and ellipticity in Table 7 
differ slightly from the values reported in Table 6. This results 
from the different way in which the Gaussian fit was applied. The 
scanning beam fit is determined by fitting the profile of Jupiter 
to timelines and limiting the fit to the data with signal-to-noise 
ratio greater than 3, while the fit of the effective beam is com¬ 
puted on GRASP maps projected in several positions of the sky 
(Planck Collaboration IV 2016). The latter are less affected by 
the noise. 


6.3. Window functions 

Window functions based on the LFI beams are needed for the 
production of the LFI likelihoods and power spectra. They are 
based on the revised FEBeCoP (effective) beams discussed ear¬ 
lier in this section, and account for the renormalization of the 
beams described in Sect. 7.2. The derivation of the 2015 window 
functions is fully described in Planck Collaboration IV (2016), 
as are the uncertainties in the window functions. The uncertain¬ 
ties are sharply reduced from the previous release and are: 0.7 % 
for the 30 GHz band (evaluated at ^ = 600); 1.0% at 44 GHz 
(also evaluated at ^ = 600); and 0.5 % in the 70 GHz window 
function at ^ = 1000. 


7.1. 4n calibration 

When we apply Eq. (1) to solve the calibration problem, we 
compute the value of 5 * D by means of a full 471 convolution 
over the sphere, between the dipole signal (plus the relativistic 
quadrupole component) and the beam response. This is different 
from what other experiments have done when using the dipole as 
a calibrator, e.g. WMAP and HFI assume the beam to be a Dirac 
delta function. Our approach allows us to properly take into ac¬ 
count the asymmetric effect of the sidelobes and the efficiency 
of the main beam during the calibration, which is critical for 
polarization, especially at low multipoles. Indeed, as discussed 
in Planck Collaboration V (2014), the introduction of 471 calibra¬ 
tion resulted in a significant improvement in the self-consistency 
of survey maps as demonstrated by null tests analysis. 

It can be demonstrated (Planck Collaboration V 2014, 2016) 
that the average level of the power spectrum, before con¬ 
volving it with the beam window function, changes with respect 
to the Dirac delta case according to the following formula: 

^ 47 : _ ~ /si “ 0sky + 0D ) 

' n )■ ® 

where /si is the sidelobe fraction of the beam, and 0/) < 0.2 %, 
0sky ~ 0.01 %, and ^ 0.01 % are quantities defined 
and discussed in Planck Collaboration V (2016); they depend 
on the beam and the scanning strategy, and they are therefore 
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Table 7. Mean and rms variation across the sky of FWHM, ellipticity, orientation, and solid angle of the FEBeCop effective beams 
computed with the GRASP beam-fitted scanning beams. Here FWHMeff is the effective FWHM estimated from the main beam solid 
angle of the effective beam, Qgff = mean(Q). 


Maps 

FWHM 

[arcmin] 

e 

•I' 

[deg] 

Q 

[arcmin' 


FWHMeff 

[arcmin] 

Frequency Maps 

LFI70GHZ .... 

13.213 + 0.034 

1.223 + 0.026 

3 + 54 

200.90 + 

0.99 

13.315 + 0.033 

LFI 44" GHz . . . 

27.000 + 0.590 

1.035 + 0.035 

0 + 50 

832.00 + 34.00 

27.100 + 0.57 

LFI30GHZ .... 

32.293 + 0.024 

1.318 + 0.037 

0 + 54 

1190.06 + 

0.69 

32.408 + 0.009 

Quadruplet Maps 

LFI 18-23 _ 

13.525 + 0.021 

1.188 + 0.021 

3 + 54 

210.13 + 

0.63 

13.618 + 0.020 

LFI 19-22 _ 

13.154 + 0.037 

1.230 + 0.027 

2 + 54 

199.19 + 

0.64 

13.259 + 0.021 

LFI 20-21 _ 

12.910 + 0.037 

1.256 + 0.036 

3 + 54 

192.58 + 

0.67 

13.037 + 0.023 

LFI 25-26 _ 

29.975 + 0.013 

1.177 + 0.030 

-2 + 47 

1019.63 + 

0.65 

29.998 + 0.009 

LFI 24 . 

23.036 + 0.014 

1.341 +0.033 

1 + 54 

603.61 + 

0.78 

23.080 + 0.015 


^ Associated errors are artificially large due to the fact that the 44 GHz maps combine the beams of the horns 24, 25, and 26, which are very 
different from each other and are located far out in the focal plane (see Fig. 2). We suggests using the quadruplet maps of LFI 25-26 and LFI 24 
separately. 


radiometer-dependent. The typical value of for the LFI 

radiometers deviates from unity by less than 1 %. 

The solution of the Eq. (1) is provided by an iterative de¬ 
striper, DaCapo, which supersedes our previous dipole-fitting 
code used in the 2013 data release. At each step the iterative 
procedure determines the radiometer gains by fitting D to the 
data, at the same time extracting the contribution from the sky 
signal. Because of the degeneracy between the overall gain level 
and the signal D, it makes sense to constrain the map dipole to 
the model. For this to work the contribution of foregrounds to 
the dipole on the sky must be included in the dipole model. 

The An beam model used in the calibration has been created 
taking into account the radiometer bandpass of each radiometer 
(measured before flight). For each radiometer about 25 realiza¬ 
tions of the main beam, intermediate beam, and sidelobe have 
been produced at fixed frequencies, chosen to fully sample the 
shape of the bandpass (as shown for the LFI 18M bandpass and 
selected frequencies in Fig. 3). Those realization were then used 
to construct a weighted An beam for each radiometer. 

7.2. Impact of 4n calibration on beam functions and source 
fluxes 

The mapping procedure assumes a pencil beam 
(Planck Collaboration VI 2016), which, in the ideal case 
of a circularly-symmetric beam, would yield a map of the 
beam-convolved sky; therefore a fraction of the signal from any 
source appears in the far sidelobes, and would be missed by 
integration of the map over the main beam alone. By the same 
token, bright resolved features in the map have temperatures 
fractionally lower than in the sky, due to signal lost in the 
sidelobes. In essence this description remains true even given 
the highly asymmetric sidelobes of the Planck beam: the main 
difference is that the far sidelobe contribution to a given pixel 
varies according to the orientation of the satellite at the time 
of observation. For LFI beams, roughly 1 % of the signal is in 
the sidelobes and this must be accounted for in any analysis of 
the maps. In particular, the flux densities of compact sources 
measured from the maps must be scaled up by the multiplicative 
factors /sour reported in Table 8. These values have been 
computed from: 

- the main beam efficiencies (Planck Collaboration IV 2016); 



iy[GHz] 

Fig. 3. Illustration of the method used to produce LFI synthetised 
beams weighted for the radiometer response (in this case LFI 
18M). The vertical lines identify the frequencies at which the 
beam has been simulated within the radiometer bandpass T(y). 
The results are then used to construct a weigthed An beam for 
each radiometer. Details on the bandpass measurements can be 
found in Villa et al. (2010). 


- a re-normalization factor introduced by the calibration 
pipeline to compensate for the missing power in the An 
beam (Planck Collaboration V 2016) (re-normalized beam 
efficiencies //norm are also reported in Table 8); 

- a factor that takes in account the horn uniform 
weights applied during the mapmaking process 
(Planck Collaboration VI 2016). 

The re-normalization factor was introduced to take in ac¬ 
count the “missing power” due to the first-order approxima¬ 
tion adopted in the computation carried out with the GRASP 
Multirefiector Geometrical Theory of Diffraction (MrGTD) soft¬ 
ware package (Planck Collaboration IV 2016). The missing 
power was proportionally distributed between main, intermedi¬ 
ate, and sidelobe parts; this procedure has an effect on the previ- 
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Table 8. Multiplicative factors that should be used to determine 
the correct flux densities from compact sources, ^renorm is the re¬ 
normalized main beam efficiency, /sour and /bi are multiplicative 
factors for flux densities and beam function (already applied in 
the delivered LFI beam functions and flux densities of sources 
in the PCCS2 (Planck Collaboration XXVI 2016)). 


^renorm /our 

Frequency Maps 


LFIVOGHz. 99.3582 1.00646 1.00346 

LFI44GHZ. 99.8827 1.00117 1.00143 

LFI30GHZ. 99.1983 1.00808 1.00258 


Quadruplet Maps 

LFI 18-23 . 99.4556 1.00547 1.00333 

LFI 19-22. 99.3764 1.00628 1.00339 

LFI 20-21 . 99.2238 1.00782 1.00368 

LFI 25-26. 99.9119 1.00088 

LFI 24. 99.8228 1.00177 


ously computed beam functions, and hence these have now been 
scaled by the factor /bi reported in Table 8. 

In practice, in order to make consistent comparison with ex¬ 
ternal data, it is essential that: 

- users interested in CMB and diffuse component analysis, 
should use the official LFI beam functions (in the LFI RIMO 
available in the Planck Legacy Archive interface^) which al¬ 
ready include the rescaling factor /bi. Alternatively, users 
who wish to perform their own beam deconvolution should 
multiply their beam functions by the factor /bi 

- users interested in point sources, the recalibration factors 
/sour should be used to obtain proper flux densities for 
sources extracted directly from LFI maps. 

7.3. Smoothing algorithm 

The uncertainty of the calibration constants increases signifi¬ 
cantly when the Planck spacecraft is aligned such that the ob¬ 
served scan circle measures a low dipole component (“minimum 
dipole”). This problem was particularly severe in the Surveys 2 
and 4, as shown in Fig. 4. 

To reduce the noise, we apply an adaptive smoothing al¬ 
gorithm that is also designed to preserve the discontinuities 
caused by abrupt changes in the working configuration of the ra¬ 
diometers (e.g. sudden temperature changes in the focal plane). 
Moreover, we apply an additive, zero-mean correction to the 
calibration constants derived from measurements of the emis¬ 
sion of an internal load kept at a stable temperature of approxi¬ 
mately 4.5 K, plus the measurement of a set of temperature sen¬ 
sors mounted on the focal plane of LFI. The amplitude of this 
correction is quite small (<^ 1 %, see Fig. 5), but its purpose is 
to account for two phenomena. 

1. During the first survey, the transponder used to download 
data to Earth was repeatedly turned on and off with a 24 h 
duty cycle. This caused periodic fluctuations in the tempera¬ 
ture of the back-end amplifiers, which were clearly traceable 
in the signal of the 4.5 K load (Mennella et al. 201 1), but are 
not visible in the calibration constants computed using the 
dipole, because of statistical noise (this is particularly true 
during dipole minima). 

^ http://archive.esac.esa.int/pla2 


2. In general, during a dipole minimum, we are not able to keep 
track of variations in the gain of the radiometers. However, 
the knowledge of the internal 4.5 K signal allows us to esti¬ 
mate an additive correction factor that mitigates the problem. 


7.4. Galactic straylight removal 

The light incident on the focal plane that does not reflect directly 
off the primary mirror (straylight) is a major source of systematic 
effects, especially when the Galactic plane intersects the direc¬ 
tion of the main spillover. This effect is now corrected by re¬ 
moving the estimated straylight signal from the timelines. To do 
this the term * Tsky of Eq. (1) has to be removed from cali¬ 
brated timelines (here ^si represents the sidelobes contribution to 
the beam). This term was computed for each radiometer by con¬ 
volving both Galactic and extragalactic emissions with the an¬ 
tenna pattern in the sidelobe region (with angle 6 > 5° from the 
main beam pointing direction). Here Tsky was estimated using 
simulated temperature and polarization maps. These included 
the main diffuse Galactic components (synchrotron, free-free, 
thermal, and anomalous dust emissions) as well as contribution 
from faint and strong radio sources and the thermal and kinetic 
Sunyaev-Zeldovich effects (although the last is barely relevant 
at LEI frequencies), as described in Planck Collaboration IX 
(2016) and Planck Collaboration X (2016). These maps are 
weighted across the band using the transmission function spe¬ 
cific to each radiometer and then summed together. Eor polariza¬ 
tion, the contributions from both synchrotron and thermal dust 
have been considered. 

The convolution was performed by transforming both the sky 
and the sidelobe pattern into spherical harmonics coefficients up 
to multipole £ = 2048. These coefficients are then properly mul¬ 
tiplied to produce an object containing convolution results for 



Fig. 4. Raw gain from radiometer 27M throughout 4 year mis¬ 
sion. Pid is a counter for pointings of the spin axis, which had 
an average duration of about 45 minutes Planck Collaboration I 
(2014). The increase of noise corresponding to the periods of 
“minimum dipole” (see text) are clearly visible for each of the 
eight surveys. Survey 2 (Pi^ range approximately 5 200-10 000) 
and Survey 4 (Pid approximately 15 700 - 20 600) exhibit a sig¬ 
nificantly higher noise, as expected from the unfavourable align¬ 
ment of the spacecraft spin axis with the Solar dipole in those 
two surveys. 
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Pid 


Fig. 5. High frequency fluctuations of the raw gain from ra¬ 
diometer 27M throughout the 4 year mission. The major de¬ 
crease in high frequency variations occurs after the transponder 
was left continuously on (at Pid = 5 000). Subsequently the high 
frequency variations are generally 1 %. 


each position on the sky (6, 0) and beam orientation angle if/. For 
each sample in the timeline, the straylight contribution has been 
evaluated by performing a polynomial interpolation. Figure 6 
shows expected Galactic straylight contribution in total inten¬ 
sity for a sample of LFI radiometers (both main and side arms), 
one at each frequency covering the full mission period. 

7.5. Colour correction 

Colour corrections are required to adjust LFI measurements for 
sources or foregrounds that do not have a thermal spectrum. Our 
initial estimates were listed in Planck Collaboration V (2014) for 
each LFI radiometer and frequency band. For power-law spec¬ 
tra they can be well approximated by a quadratic relationship 
between flux density and spectral index a (or equivalently tem¬ 
perature spectral index p = tr - 2), where the quadratic coeffi¬ 
cient is proportional to the square of the fractional bandwidth, 
and the linear term mainly depends on the value of the chosen 
reference frequency (Leahy & Foley 2006). The constant com¬ 
ponent is constrained by the requirement of zero colour correc¬ 
tion for the CMB spectrum, so there are two free parameters in 
the model. Accurate quadratic fits are used in the fastcc IDL 
code included in the Planck unit conversion and colour correc¬ 
tion software package. 

The more detailed component separation analysis for the 
2015 release (Planck Collaboration X 2016) has allowed us to 
further constrain the colour corrections, which in the 2013 re¬ 
lease were based purely on ground-based measurement and 
modelling of the radiometer bandpasses. In recent analyses, we 
used separate maps from each of the three co-scanning pairs 
of 70 GHz horns. The analysis uses maps from LFI, HFI, and 
WMAP, which includes several pairs of channels spaced closely 
in frequency. Using the nominal colour corrections for the three 
instruments, highly significant and systematic residuals were 
found to our best-fit models for the strong Galactic emission, 
which resemble gain errors; however, gain errors can be ruled 
out, because there were no detectable residuals correlated with 


Table 9. Coefficients for parabolic fits to the LFI colour correc¬ 
tions C(a), revised from the 2013 values, based on the band¬ 
pass shifts derived by Commander component separation code 
(Planck Collaboration X 2016). 


Horns Band vq cq Ci C 2 

[GHz] [GHz] 

27,28 . 30 28.4 1.005 0.0030 -0.0030 

24,25,26... 44 44.1 0.995 0.0060 -0.0017 

18,23 . 70 70.4 0.983 0.0142 -0.0032 

19,22 . 70 70.4 1.010 -0.0007 -0.0033 

20,21 . 70 70.4 0.977 0.0176 -0.0031 

All. 70 70.4 0.990 0.0107 -0.0033 


the CMB emission. We thus assume that the previous colour cor¬ 
rections caused the residuals, and have tried to improve them. 

A first attempt has been made to derive improved colour cor¬ 
rections by fitting for a frequency shift in the bandpass as part 
of the component separation analysis. This minimal model was 
adopted to avoid a strong degeneracy between the bandpass re¬ 
calibration and the foreground spectral models; it is certainly an 
oversimplification. The resulting fractional change of frequency 
is 1.0 + 0.3 % at 30 GHz, 0.2 + 0.2% at 44 GHz, and -0.6%, 
1.6% and 0.7% (all +1.4%) for the three 70 GHz horn pairs 
(18 and 23, 19 and 22, 20 and 21, respectively). The uncertain¬ 
ties quoted here are the absolute ones. For convenience. Table 9 
lists the parameters of our parabolic fit to the colour corrections 
derived from the shifted bandpasses, where for a map thermody¬ 
namic temperature f, the Rayleigh-Jeans brightness temperature 
at the reference frequency vq is given by 

T(yo)[KRj] = f [Kcmb] ^Ar(vo) C(a), ( 3 ) 

where ^Ar(y) = the coefficients in Table 9 

give the colour correction as C(a) = cq + cia + C 2 cr^. Because 
they are based on a simplified analysis, these values should be 
treated with some caution; the revised colour corrections have 
only been tested for spectral indices near that of the dominant 
foregrounds at each frequency, namely -1 a ^ 0at30 and 
44 GHz, and 0 ^ a ^ 2.5 at 70 GHz. We plot the old and new 
corrections in Fig. 7. 

In Planck Collaboration V (2014) we gave a rough indirect 
estimate of the expected uncertainties in the colour corrections, 
assuming that the errors for the individual radiometer band- 
passes were uncorrelated. The revised corrections at 30 GHz dif¬ 
fer from our original ones by 2-3 % for for a -1, and this 
change is almost an order of magnitude larger than our origi¬ 
nal error estimate. In retrospect, our assumption of uncorrelated 
errors was flawed for this particular channel, since our ground- 
based estimates of the bandpass shape were particularly sensitive 
to modelling assumptions. This arose because the bands still had 
significant response at the low-frequency end of the the directly- 
measured range. As a result, it seems likely that the actual differ¬ 
ence between our 2013 model and the true 30 GHz bandpasses 
is more in the nature of an upward revision of the low-frequency 
cutoff than a uniform shift to higher frequency; if this is the case, 
then our 2015 estimate for C(a) will still show too much curva¬ 
ture in this band. 

The estimated bandpass shifts at 44 and 70 GHz are not 
significant, but they correspond to colour-correction changes 
of 0.5 % and 1 %, respectively, consistent with our original er¬ 
ror estimates; moreover, as explained in Planck Collaboration X 
(2016), the relative shifts between the 70 GHz horn pairs are 
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Fig. 6. Simulated Galactic stray light in total intensity for representative LFI radiometers for the full mission period. Top: 70 GHz 
radiometer 18M (right) and 18S (left). Middle: 44 GHz radiometer 24M (right) and 24S (left). Bottom: 30 GHz radiometer 27M and 
27S (left). The faint stripes paralleling the scanning direction are due to the different coverage of the sky during different surveys. 



a 


Fig. 7. Colour corrections C(a) versus intensity spectral index 
a. Solid lines are the current corrections given by Table 9, while 
dashed lines are the 2013 values. Red curves are for the 30 GHz 
band, green for 44 GHz, and blue for 70 GHz. Note that the cor¬ 
rections have only been validated for a ^ 0 at 30 and 44 GHz, 
and for 0 < a < 2.5 at 70 GHz. 


known much more accurately than their absolute values, and 
are certainly important. We therefore recommend the use of the 
revised colour corrections listed here. The uncertainties in the 
correction should be taken to be approximately |yS| x 0.3 % for 
ah channels, as long as the spectral index is close to the well- 


sampled range, -3 ^ ^ 1. We note that, by construction, the 

colour correction tends to unity for emission having the colour 
of the CMB, and so C remains accurately equal to unity when 
P = PcMR • 

Since it is possible to make total intensity sky maps from 
the data for each individual LFI feed horn (averaging the data 
from the M and S radiometers), it will be possible to improve 
the colour corrections individually for each horn, and we plan to 
do that for the next release. 

7.6. Summary of changes in LFI calibration 

In this subsection, we summarize the changes in the overall cal¬ 
ibration of the PlanckLFl channels that have resulted from dif¬ 
ferent procedures adopted since Planck Collaboration II (2014) 
and from our deeper understanding of instrumental systematics 
and their effect on calibration. 

- Overall calibration. Improved accounting for beam effects 
and other changes discussed in Sect. 6 and Sect. 7 produces 
a small upward shift in the calibration for the three LFI chan¬ 
nels. In addition, our current use of the orbital dipole for 
the determination of the solar dipole used for calibration has 
shown that the previous calibration based on the WMAP so¬ 
lar dipole was 0.28 % low for ah frequencies. Combining 
these effects, we hnd the following upward shifts in LFI cal¬ 
ibration: 0.83%, 0.72%, and 0.95% for 30 GHz, 44 GHz, 
and 70 GHz, respectively. 

- Uncertainties in calibration. Improved understanding and as¬ 
sessment of the impact of various systematic effects on cal- 


11 





Planck Collaboration: LFI data processing 


ibration have allowed us to refine our estimates of overall 
calibration uncertainty. The uncertainties are 0.35 %, 0.26 %, 
and 0.20 % for 30, 44, and 70 GHz, respectively. 

- Window function. We now use 4?: beams, rather than a pen¬ 
cil beam approximation. LFI window functions properly take 
account of the small amount of missing power in the side- 
lobes (a roughly 0.4 % effect at most, see Table. 8). 

“ Flux densities of compact sources. Our current use of a 471 
beams also means that flux densities of compact sources 
need to be boosted by a small factor if they are derived 
from the LFI maps (again, see Table 8). Flux densities in 
the PCCS2, on the other hand, are already corrected for this 
factor. 

8. LFI noise estimation 

8.1. Radiometer noise model 

A detailed knowledge of instrumental noise properties is funda¬ 
mental for several stages of the data analysis. First of all evolu¬ 
tion in time of basic noise properties (e.g. white noise variance) 
throughout the entire mission lifetime is an important and sim¬ 
ple way to track possible variations and even anomalies in the 
instrument behaviour. In addition, noise properties serve as in¬ 
puts for the Monte Carlo noise simulations (used, e.g. for power 
spectrum estimation) and also give correct weights for properly 
combining different detectors. 

We proceed as already shown in Planck Collaboration II 
(2014) using an implementation of a Monte Carlo Markov chain 
(MCMC) approach to estimate basic noise properties. As before, 
the noise model is 


P(f) = 


1 -r 



..^ee / 


(4) 


where cr^ is the white noise level, and /knee and p describe the 
non-white component of the instrumental noise. To evaluate cr^, 
we take the mean of the noise spectrum in the last few (typically 
10 %) of the bins at the highest frequency, which exhibits a fiat, 
high-frequency tail, as shown in Figure 8. At 30 GHz the knee- 
frequency is /knee ~ 100 mHz and therefore a smaller percentage 
of data has been taken for computing cr^. These values for a are 
given in Table 10. Once this is done we can proceed with the 
evaluation of the other two parameters. After discarding a burn- 
in period from our chains, we obtained the best-fit and variances 
values reported in Table 11 . 


8.2. Updated noise properties 

We estimate noise properties at the radiometer level using the 
MCMC approach. As already done with the previous data re¬ 
lease, we work with calibrated data and select chunks of data 
5 days long and process them with the roma generalized least- 
squares mapmaking algorithm (de Gasperis et al. 2005). The 
outputs are frequency spectra that are then fitted for the basic 
noise parameters. Results are summarized in Tables 10 and 11, 
for the white noise sensitivity and 1 // noise parameters, respec¬ 
tively. These numbers are the medians, computed from the fit 
results throughout the whole mission lifetime. 

Time variations of the noise properties are a good indicator 
of possible changes in instrument behaviour. There are known 
events that caused such variations, such as the sorption cooler 
switchover at OD 460 (Planck Collaboration I 2014). Indeed, 
variations in noise properties due to changes in temperature are 


Table 10. White noise sensitivities for the LFI radiometers. 



White Noise Sensitivity 


Radiometer M 
[juKcmb 

Radiometer S 
[yuKcMB s‘/2] 

70 GHz 

LFI 18. 

513.0+2.1 

467.2+2.3 

LFI 19. 

579.6+2.2 

555.0+2.2 

LFI 20. 

587.3+2.1 

620.5 +2.7 

LFI 21. 

451.0+1.7 

560.1+2.0 

LFI 22. 

490.8+1.5 

531.3+2.3 

LFI 23. 

504.3+1.8 

539.7+1.8 

44 GHz 

LFI 24. 

463.0+1.4 

400.7+1.3 

LFI 25. 

415.3+1.5 

395.4+2.9 

LFI 26. 

483.0+1.9 

423.2+2.5 

30 GHz 

LFI 27. 

281.5+2.1 

303.2 + 1.8 

LFI 28. 

317.7+2.4 

286.5 +2.3 


expected as the performance of the first cooler degraded, as well 
as when the second cooler came in and took time to stabilize 
the temperature. Figure 8 shows a sample of noise spectra for 
radiometers LFI27M, LFI25S, and LFI18S, spanning the whole 
mission lifetime. The white noise is stable at the level of 0.3%. 
As already noted in the previous release, knee-frequencies and 
slopes are stable until OD 326 and show significant variations af¬ 
terwards, altering the simple “one slope, one knee” model. This 
is due to the progressive degradation of the first sorption cooler 
and the insertion of the second one. Once the environment be¬ 
came thermally stable, the spectra moved back towards their ini¬ 
tial shape. Of course this is evident at different levels in the in¬ 
dividual radiometers, depending on their frequency, position on 
the focal-plane, and susceptibility to thermal instabilities. 


9. Mapmaking 

Mapmaking is the last step in the LFI pipeline, after calibration 
and dipole removal, and before bandpass correction and com¬ 
ponent separation. Mapmaking takes as its input the calibrated 
timelines, from which the 47r convolved dipole and Galactic 
straylight signal has been removed. Output consists of sky maps 
of temperature and Q and U polarization, and a description of 
the residual noise in them. 

An important part of the mapmaking step is the removal of 
correlated 1// noise. An optimal mapmaking method will re¬ 
move the noise as accurately as possible, while simultaneously 
keeping systematics at an acceptable level. 

LFI maps were produced by the Madam mapmaking code 
Keihanen et al. (2005). The code is the same as used in the 2013 
release. In the following we give a short overview, and point 
out aspects relevant to polarization (see Planck Collaboration VI 
(2016) for details). 

Madam removes the correlated noise using a destriping tech¬ 
nique. A noise prior is used to improve the map quality further. 
The correlated noise component is modelled by a sequence of 
baseline offsets. The choice of the baseline length is a trade-off 
between computational burden and optimal noise removal. We 
have chosen to use 1 s long baselines for 44 and 70 GHz, and 
0.25 s for 30 GHz where the typical knee frequencies are higher. 
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Table 11. Knee frequencies and slopes for the LFI radiometers. 



Knee Frequency /knee [mHz] 

Slope p 


Radiometer M 

Radiometer S 

Radiometer M 

Radiometer S 

70 GHz 

LFI 18. 

14.8 + 2.5 

17.8 + 1.5 

-1.06 + 0.10 

-1.18 + 0.13 

LFI 19. 

11.7 + 1.2 

13.7 + 1.3 

-1.21+0.26 

-1.11+0.14 

LFI 20. 

8.0+1.9 

5.7+ 1.5 

-1.20 + 0.36 

-1.30 + 0.41 

LFI 21. 

37.9 + 5.2 

13.3 + 1.5 

-1.25 + 0.09 

-1.21+0.09 

LFI 22. 

9.7+ 2.3 

14.8 + 6.7 

-1.42 + 0.23 

-1.24 + 0.30 

LFI 23. 

29.7 + 1.1 

59.0+1.4 

-1.07 + 0.03 

-1.21+0.02 

44 GHz 

LFI 24. 

26.8 + 1.3 

88.3 + 8.9 

-0.94 + 0.01 

-0.91 + 0.01 

LFI 25. 

20.1+0.7 

46.4 + 1.8 

-0.85 + 0.01 

-0.90 + 0.01 

LFI 26. 

64.4+1.9 

68.2 + 9.5 

-0.92 + 0.01 

-0.76 + 0.07 

30 GHz 

LFI 27. 

174.5 + 2.9 

108.8 + 2.5 

-0.93 + 0.01 

-0.91+0.01 

LFI 28. 

130.1+4.4 

43.1 + 2.4 

-0.93 + 0.01 

-0.90 + 0.02 




Fig. 8. Noise spectra throughout the mission lifetime for a 70 GHz radiometer 18M (left), 25S (44 GHz; middle), and 27M (30 GHz; 
right). Spectra are shown for the ranges from OD 100 (blue) to OD 1526 (red), spaced about 20 ODs apart. White noise is stable 
at the level of 0.3%, while low-frequency noise shows variations both in slope and knee-frequency, with different amplitude for 
different radiometers. 


The full time-ordered data stream is modelled as 

y = Pm -\-Pa-\- n. (5) 

Here vector a represents the baselines, and F is formally a ma¬ 
trix that spreads the baselines into time-ordered data. Vector n 
represents white noise, and P is a pointing matrix that picks a 
time-ordered data stream from the sky map m. Map m has three 
columns, corresponding to the three Stokes components 1, Q, 
and U. 

The noise prior describes the expected correlation between 
baseline amplitudes, 

c„ = (aa'T). (6) 

The prior is constructed from the known noise parameters pre¬ 
sented in the previous section (knee frequency, white noise 
sigma, and spectral slope). The noise prior provides an extra 
constraint which makes it possible to extend the destriping tech¬ 
nique to very short baseline lengths, allowing for more accurate 
removal of noise. 

With the assumptions above, the baseline vector a can be 
solved from the linear system of equations 

(FTc;'ZF + c-')a = FTc;‘Zy, (7) 

where 

z = i-P(pTc;‘Pr‘PX‘- (8) 


Here Cw is a diagonal weighting matrix. The final map is con¬ 
structed as 

m = (pTc;‘Pr‘P'fC;‘(>- - Fa). (9) 

The destriping technique constructs the final map through 
a procedure in which one first solves for the baselines, and then 
bins the map from the data stream from which the baselines have 
been removed. This two-step procedure provides a way of re¬ 
ducing systematics. We control the “signal error” by applying a 
mask in the destriping phase, while still binning the final map to 
cover the whole sky. Signal error is the uncertainty in baseline 
determination that arises from deviations of the actual sky signal 
from the model Pm. The main sources of signal error are signal 
variations within a pixel, differences in radiometer frequency re¬ 
sponses (bandpass mismatch), and beam shape mismatch. The 
error arises mainly at low Galactic latitudes, where signal gradi¬ 
ents are strong. 

The choice of the destriping mask is a trade-off between ac¬ 
ceptable signal error level and noise removal. A mask that is too 
wide may lead to a situation where there are not enough cross¬ 
ing points between scanning rings to reliably determine the noise 
baselines. 

It can be shown that residual noise is minimized when Cw 
equals the variance of white noise in time-ordered data. In or¬ 
der to reduce leakage from temperature to polarization, however. 
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we apply horn-uniform weighting, which differs from this ideal 
case. We replace the white noise variance by the average of the 
variances of the two radiometers of the same horn. This has the 
effect that the systematic error related to beam shape mismatch, 
which is strongly correlated between the radiometers, largely 
cancels out in polarization analysis. Thus we are reducing the 
leakage from temperature to polarization. 

Along with maps of the sky. Madam provides a covariance 
matrix for residual white noise in the maps. This consists of a 
3x3 matrix for each pixel, describing the correlations between 
/, 2, and U components in the pixel. White noise is uncorre¬ 
lated between pixels. Correlated noise residuals are captured by 
the low-resolution noise covariance matrix describe in Sect. 9.3 
below. 

Madam produces its output maps in HEALPix format 
(Gorski et al. 2005). For the bulk of the products we used res¬ 
olution Aside = 1024, and the same resolution was used when 
solving the destriping equation. Maps at 70 GHz were also pro¬ 
duced with Aside = 2048. 

To accurately decompose the map into /, 2, and U compo¬ 
nents it is necessary to have several measurements from the same 
sky pixel, with different parallactic angles. If this is not the case, 
the pixel in question is eliminated from analysis. Madam uses as 
rejection criterion the reciprocal condition number of the matrix 
P^C;ip. 

9.1. Low-resolution data set 

Low-resolution products are an integral part of the low-L like¬ 
lihood. To fully exploit the information contained in the largest 
structures of the microwave sky, a full statistical description of 
the residual noise present in the maps is required. This infor¬ 
mation is provided in the form of pixel-pixel noise covariance 
matrices (NCVMs). However, due to resource limitations they 
are impossible to employ at native map resolution. Therefore a 
low-resolution data set is needed for the low-^ analysis; this data 
set consists of low-resolution maps and corresponding noise co- 
variance matrices. At present, the low-resolution data set can be 
efficiently used only at resolution Aside = 16, or lower. All the 
low-resolution products are produced at this target resolution. 

9.2. Low-resolution maps 

The low-resolution maps, shown in Figure 9, are constructed by 
downgrading the high-resolution maps (described in the previ¬ 
ous section) to the target resolution. We chose to downgrade the 
maps using a “noise-weighted” scheme. 

The noise-weighted scheme has also been used in previ¬ 
ous studies (see, e.g. Planck Collaboration II 2014). The noise- 
weighted map corresponds to a map that is first destriped at the 
high resolution, and the destriped TOI is directly binned onto 
the low target resolution. This approach gives adequate con¬ 
trol over signal and noise in the resulting maps. However, con¬ 
cerns have been raised that the noise-weighted scheme transfers 
signal from one pixel to another. As a consequence we em¬ 
ploy Gaussian smoothing to minimize this effect, at the cost 
of some increase in noise. After downgrading, the temperature 
component is smoothed with a Gaussian window function with 
FWHM = 440'. We will re-examine this choice in the next re¬ 
lease. 

In practice the high-resolution maps are noise-weighted to 
an intermediate resolution of Aside = 32. The Stokes I part of the 
map is expanded in spherical harmonics, the expansion is treated 


with the smoothing beam, and the final map is then synthesized 
at the target resolution. The last step of resolution downgrading 
for Stokes Q and U maps, however, is performed by carrying out 
naive averaging of higher resolution pixels. 

Due to the chosen downgrading scheme the resulting NCVM 
will be singular. We regularize the problem by adding some 
white noise both to the maps and matrices. Specifically we add 
2 jiiK for /, and 0.02 juK for Q and U at Aside = 16 resolution. 

9.3. Noise covariance matrices 

The statistical description of the residual noise present in a low- 
resolution map is given in the form of a pixel-pixel noise covari¬ 
ance matrix, as described in Keskitalo et al. (2010). The NCVM 
formalism describes the noise correlations of a map produced at 
the same resolution as the noise covariance matrix. Therefore, 
for an exact description we should construct the matrices at res¬ 
olution Aside = 1024 and subsequently downgrade to the target 
resolution. This is computationally impractical. Therefore the 
matrices are computed at the highest possible initial resolution, 
and then downgraded to the target resolution. For consistency the 
noise covariance matrices must go through the same processing 
steps as applied to the low-resolution maps. 

The Madam/TOAST code, a Time Ordered Astrophysics 
Scalable Tools (TOAST) port of Madam , was used to produce 
the pixel-pixel noise covariance matrices (Keihanen et al. (2010) 
and http://tskisner.github.io/TOAST/. The TOAST in¬ 
terface was chosen on the basis of added flexibility and speed; 
see (Planck Collaboration VI 2016). 

The outputs of Madam/TOAST software are inverse-noise co- 
variance matrices, specifically one inverse matrix per radiome¬ 
ter for a given time period. Because inverse NCVMs are addi¬ 
tive, the individual inverse matrices are merged together to form 
the actual inverse NCVM. To obtain the noise covariance ma¬ 
trix from its inverse, the matrices are inverted using the eigen- 
decomposition of a matrix. These intermediate-resolution matri¬ 
ces are then downgraded using the same downgrading scheme 
as applied to the maps. The matrices are regularized by adding 
the same level of white noise to the diagonal elements of the 
covariance matrix as to the low-resolution maps. 

The noise covariance matrix computation takes two inputs: 
the detector pointing; and noise estimates. Since the matrices are 
calculated with Madam/TOAST, we use the pointing solution pro¬ 
vided by TOAST. For more details see Planck Collaboration XII 
(2016). We also use the most representative noise model avail¬ 
able, namely the FFP8 (full focal plane 8 simulations) noise es¬ 
timates (Planck Collaboration XII 2016). The noise model com¬ 
prises daily 1// model parameters. 

The key parameter in the NCVM production is the 
baseline length. We have demonstrated in an earlier study 
(Planck Collaboration II 2014) that using shorter baseline 
lengths when producing the noise covariance matrix produc¬ 
tion better models the residual noise. Therefore we chose to use 
0.25 s baselines for the 30 GHz LFI frequency channel; we show 
in (Planck Collaboration VI 2016) that 1.0 s is adequate for the 
44 GHz and 70 GHz channels. Reducing the baseline length still 
further gives only a marginal improvement, while the resource 
requirements increase rapidly. 

Previous studies (Planck Collaboration II 2014) have also 
shown that matrices should be calculated at the highest com¬ 
putationally feasible resolution. For the current release the ini¬ 
tial resolution is Aside = 64. Increasing the initial resolution be¬ 
yond Aside = 64 is likely to improve results, but the matrix size 
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Fig. 9. LFI full mission low-resolution maps, A^side = 16. From left to right 30 GHz, 44 GHz, and then 70 GHz: top intensity /; 
middle polarization Q component; and bottom polarization U component. Units are pY^cMB • 





will be 16 times larger, i.e. 2.5 TB. Inverting such a matrix is a 
formidable task. 

The noise covariance computation makes two further devia¬ 
tions from the high-resolution mapmaking: it does not take into 
account the destriping mask; and the horns are not uniformly 
weighted. The effect of these differences is much smaller than is 
obtained by either decreasing the baseline length or increasing 
the destriping resolution in the production. For more details see 
Planck Collaboration VI (2016). 


10. Overview of LFI map properties 

Figures 10 to 12 show the 30, 44, and 70 GHz frequency maps 
created from LFI data. The top panel in each figure is the tem¬ 
perature (/) map, based on the full observation period at native 
resolution and HEALPix Vside = 1024. The middle panel is the 
Q polarization component, while the bottom panel is the U po¬ 
larization component at Vside = 256 smoothed at 1° resolution. 
In Fig. 13 the eight surveys at 30, 44, and 70 GHz are shown; 
the grey areas identify the regions of the sky not observed in 
each survey. Table 12 reports the main parameters used in the 
mapmaking process. 

The delivered maps have been processed in order to remove 
any spurious zero-level (or monopole term). To do this we im¬ 
plemented the following procedure. We derived from LFI data 
only an estimation of the CMB signal by processing 1° smoothed 
maps with an ILC (Internal Linear Combination) method, as 
described in Eriksen et al. (2004). We then smoothed the sin¬ 
gle frequency LFI maps at the same resolution and subtracted 
the CMB estimate. For each map we used the variation with 
Galactic latitude of the remaining Galactic emission signal to es¬ 
timate the zero-level. We assumed a simple plane-parallel model 
for Galaxy emission and fit the data with a functional form as 
T = Acsch B in the range -90° < h < -15°, using the same 
mask as employed in the mapmaking procedure. The value of B 
is the zero-level we are looking for, which has to be subtracted 


from the maps in order to obtain an overall “null” zero-level. 
This value is reported in Table 12. 

Finally Table 13 lists the delivered maps along with the data 
period used to create them. All have HEALPix resolution Aside = 
1024; in the case of 70 GHz we also provide maps at a higher 
resolution. Aside = 2048. 


11. Polarization 

The most important new results in this release are polarization 
measurements. The maps of Stokes Q and U at each LFI fre¬ 
quency are shown in Figs. 10, 11, and 12 at 30, 44, and 70 GHz, 
respectively. The 70 GHz polarized data play a critical role in 
the construction of the Planck low-^ likelihood, as described 
in Planck Collaboration XI (2016). Given the small amplitude 
of CMB polarization, we have paid careful attention to system¬ 
atic effects that could bias our polarization results. The dom¬ 
inant effect is leakage of unpolarized emission into polariza¬ 
tion (Leahy et al. 2010), which we describe in detail in section 
11.2. An overview of systematics impacting both temperature 
and polarization data is provided in section 12.6, while a full 
account of the 2015 LFI systematic error budget is given in 
Planck Collaboration III (2016). 

11.1. Bandpass mismatch leakage 

Any difference in gain between the two arms of an LFI radiome¬ 
ter will result in leakage of unpolarized emission into the po¬ 
larization signal. Since gains are calibrated by observations of 
the CMB dipole, exact gain calibration would ensure that un¬ 
polarized, well resolved, CMB emission perfectly cancels in the 
polarization signal. 

However, because the bandpasses of the two arms are not 
identical, unpolarized foreground emission, if it has a different 
spectrum from the CMB, will still appear with different ampli¬ 
tudes in the two arms and therefore leak into polarization. This is 
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Fig. 10. LFI maps at 30 GHz. Top: intensity 1. Middle: polarization Q component. Bottom: polarization U component. Polarization 
components are at Nside = 256 and smoothed at 1°, the intensity is left at the native Nside = 1024. Units are yuKcMB- The polarization 
components have been corrected for the bandpass leakage effect (see Sect. 11). 
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Fig. 11. Same as Fig. 10 for 44 GHz. 
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Fig. 12. Same as Fig. 10 for 70 GHz. 
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Fig. 13. Individual survey temperature maps. Left: 30 GHz. Middle: 44 GHz. Right: 70 GHz. From top to bottom are Surveys 1 to 8. 
the gray area identify the regions of the sky not observed in each survey that depends from the spin axis orientation. 
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Table 12. Frequency-specific mapmaking parameters and related information. Details are reported in Planck Collaboration VI 
(2016). 


Channel 

/samp [Hz]^ 

Baseline length*^ 

Resolution^ 

Monopole, B^ 

[S] 

Samples 

Vide 

[arcmin] 

[//Kcmb] 

30 GHz 

32.508 

0.246 

8 

1024 

3.44 

-83.5 + 0.7 

44 GHz 

46.545 

0.988 

46 

1024 

3.44 

-40.6 + 0.7 

70 GHz 

78.769 

1.000 

79 

1024/2048 

3.44/1.72 

-35.7 + 0.6 


^ Sampling frequency. 

^ Baseline length in seconds and as a number of samples. 

^ HEALPix Vside resolution and pixel averaged size. 

^ Monopole removed from the maps; the value is included in the header fits. 


Table 13. Periods covered by the released maps. 


Period 

OD range^ 

30 GHz 


44 GHz 


70 GHz 

Sky cov. [%] 

Horns’^ 

Sky cov. [%] 

Horns^ 

Sky cov. [%] 

Homs'' 

Full. 

[ 91-1543] 

100.00 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Survey 1. 

[91- 270] 

97.20 

27,28 

93.93 

24,25,26 

97.94 

18,19,20,21,22,23 

Survey 2. 

[ 270- 456] 

97.48 

27,28 

93.31 

24,25,26 

97.47 

18,19,20,21,22,23 

Survey 3. 

[ 456- 636] 

97.62 

27,28 

93.65 

24,25,26 

97.61 

18,19,20,21,22,23 

Survey 4. 

[ 636- 807] 

91.88 

27,28 

89.53 

24,25,26 

92.40 

18,19,20,21,22,23 

Survey 5. 

[ 807- 993] 

90.89 

27,28 

88.43 

24,25,26 

92.44 

18,19,20,21,22,23 

Survey 6. 

[993-1177] 

87.79 

27,28 

86.10 

24,25,26 

89.95 

18,19,20,21,22,23 

Survey 7. 

[1177-1358] 

85.40 

27,28 

83.70 

24,25,26 

88.43 

18,19,20,21,22,23 

Survey 8. 

[1358-1543] 

80.01 

27,28 

78.92 

24,25,26 

83.83 

18,19,20,21,22,23 

Special survey^ . . . . 

. . Sl+S3+S[5-8] 

100.00 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Year 1. 

[ 91- 456] 

100.00 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Year 2. 

[ 456- 807] 

99.98 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Year 3. 

[807-1177] 

99.72 

27,28 

99.91 

24,25,26 

99.65 

18,19,20,21,22,23 

Year 4. 

[1177-1543] 

95.67 

27,28 

96.87 

24,25,26 

97.38 

18,19,20,21,22,23 

Year 1 - Year 2^ ... . 

Y1+Y2 

100.00 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Year 1 - Year 3^ . . . . 

Y1+Y3 

100.00 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Year 2 - Year 4^ . . . . 

Y2+Y4 

100.00 

27,28 

100.00 

24,25,26 

100.00 

18,19,20,21,22,23 

Year 3 - Year 4^ . . . . 

Y3+Y4 

99.77 

27,28 

99.90 

24,25,26 

99.66 

18,19,20,21,22,23 


^ OD (operational day) is defined as the time period between one daily telecommand and the succeeding one; it corresponds to about 24 h. 

^ Full period maps have been delivered for frequency, pairs of horn (18-23, 19-22, 20-21, 25-26), single radiometers, and in the case of horn 24, 
for a single horn. 

^ This special survey has been created excluding Surveys 2 and 4 for low multipole analysis, see Sect. 13. 

^ Year differences. 


“bandpass mismatch” leakage, which was discussed extensively 
in Leahy et al. (2010). 

In principle, two approaches can be used to correct for it. The 
first exploits the fact that the bandpass leakage is independent of 
the polarizer orientation, and performs a “blind” separation us¬ 
ing observations of a given pixel with multiple orientations of the 
same radiometer. With the second method, we can characterize 
both the instrumental bandpass mismatch, and the foreground 
spectrum and intensity, and hence predict the leakage explicitly. 

The blind approach was used by WMAP (Page et al. 2007), 
but for most sky pixels it is not effective for Planck, because only 
a relatively small range of detector orientations are available; this 
causes very large covariances between the leakage and the true 
Q and U values, effectively increasing the Q and U noise by a 
large factor. Hence we use the predictive method to calculate the 
leakage in our Q and U maps, and subtract it. We discuss in turn 
the determination of the foreground model, the derivation of the 
instrumental term, and the algorithm for making the correction. 


11.2. Leakage maps 

The spectra of all important LFI foregrounds are very smooth 
continua, and so to a good approximation can be modelled as a 
power law within the bandpass at each LFI frequency band. As 
described by Leahy et al. (2010), the leakage into the polariza¬ 
tion signal recorded by radiometer k can be written as 

Sk = ClkiP - PcMB)Ty^ , (10) 

where the tz-factor characterizes the bandpass mismatch (see 
next subsection), (3 = d\nT^ jd In y\y^ is the spectral index of the 
foreground within the band and is the foreground Rayleigh- 
Jeans brightness temperature at the band fiducial frequency vq. 
We separate this into the instrumental ( 2 -factor, and an astro- 
physical leakage term L = {f3 -yScMB)^^^- We derive L from 
our Bayesian component separation analysis, as described in 
Planck Collaboration X (2016). The analysis incorporates the 
Planck full-mission data, along with the WMAP 9-year maps 
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and the Haslam et al. (1982) 408 MHz map, to give 15 data 
points at each pixel. This was an earlier run than the one de¬ 
scribed in Planck Collaboration X (2016): the Planck maps were 
from a slightly earlier version of the calibration pipeline; the 
original bandpass models were used to make colour correc¬ 
tions; only a single spinning dust component was included in the 
model, not two; and the synchrotron template from the Galprop 
code was scaled only in amplitude, not in frequency. 

This analysis produces numerous Gibbs-sampled realiza¬ 
tions of the astrophysical component parameters, from which 
and p can be reconstructed at any given frequency, for each pixel 
in each realization, j. In practice, we evaluate these individually 
for each component /, to find a leakage map for each compo¬ 
nent, U., and then sum the components to give Lj = 2/ kJj. This 
is not only more straightforward to evaluate, but also automat¬ 
ically corrects for any in-band spectral curvature caused by the 
superposition of foregrounds with similar amplitudes but very 
different spectral indices. The final leakage map is then simply 
the average over the realizations, L = {IIN) Lj. In practice 
we use 1000 realizations taken after the sampling chains have 
successfully burnt in. 

The uncertainty in L at each pixel is based on the variance 
over the Gibbs realizations, ctl. However, we also have a mea¬ 
sure of goodness-of-fit of the model measured per pixel for 
each realization. Because our MCMC chains are well burned in, 
we work with the average = (^/N) YjX] over all realizations 
(but still separate for each pixel). In regions of strong foreground 
emission the component separation is limited not by noise but by 
a mismatch between the assumed algorithmic form of the model 
and the actual spectrum, signalled by high;^^. Because the model 
is non-linear and many of the model parameters are subject to 
strong prior constraints, the;^^ statistic is not expected to follow 
a distribution with the number of degrees of freedom equal 
to the number of data points. We therefore define a fiducial x^ 
equal to the median x^ over the whole sky, which of course is 
dominated by the high-latitude regions, where the foregrounds 
are weak, and therefore the component separation residuals are 
dominated by noise. We adopt an empirical correction to the un¬ 
certainty by multiplying ctl by the square root of the ratio of the 
mean;^^ to our fiducial value wherever this ratio exceeds unity. 

The component separation analysis must be done at identi¬ 
cal resolution for all frequency channels, and this was chosen as 
1° FWHM to allow use of the 408-MHz survey. Consequently, 
polarization maps corrected for bandpass mismatch leakage 
are only available at this or lower resolution. Since the full- 
resolution polarization maps have a signal-to-noise ratio of much 
less than unity for nearly all pixels, most scientific analysis must 
in any case be done with smoothed maps, or equivalently with 
only the low multipoles in harmonic space, so the low resolu¬ 
tion of the leakage maps is not a problem for most purposes. 
However, full resolution data are needed to give the most accu¬ 
rate polarimetry of point sources. A special procedure was there¬ 
fore used to correct the polarization of sources, as described in 
Planck Collaboration XXVI (2016). 


11.3. a-f actors 

Ground-based measurements of the LFI instrumental band- 
passes (Zonca et al. 2009) are not accurate enough for our pur¬ 
pose. Fortunately, to a good approximation, the bandpass mis¬ 
match can be characterized by a single parameter, the a-factor, 
which quantifies the difference in effective frequency between 


the two bandpasses: 


VefF,s Vgff jji 

2yo 


(11) 


where “s” and “m” refer to the side and main arms of the ra¬ 
diometer, respectively, and vq = (veff,s + yefF,m)/2. We deter¬ 
mine the ^-factors from flight data by first using the blind ap¬ 
proach at each frequency to estimate (I, Q,U,S 1,82 - - •) (here¬ 
after IQUSS) Sit each pixel, where Skis the spurious signal from 
each RCA. Taking the 30 GHz data as an example, the S k maps 
are defined as 


( dsi= IQcos(2il/si) + Usiniliksi) + , 

\dmi= I 2cos(2j/rjni) -r Usm(2i//mi) -Si, 

T pT 98 / ^ + QcOS(2i//s2 ) + U Sin(2j/rs2) + ‘S' 2 , 

\dm2=I^ QC0S(2lkm2) + Usm(2lkm2) " ^2 • ^ ^ 

This can also be written in a more compact form 

df = I Qcos(2if/i) + Usixi(2if/i) + triiS i -i- 0^28 2 ? (15) 


where ai and a 2 take the values -1,0, -rl, depending on the ra¬ 
diometer. To estimate m = [I, Q, U, 81 , 82 ] we need to solve a 
problem similar to mapmaking, where the noise covariance ma¬ 
trix per pixel is given by the usual 3x3 matrix block from 
Madam, with two more columns (and rows) in the form 

Mp = X 

i€p 


'... ... ... ai a2 

ai cos(2i//i) a 2 cos(2i//i) 
ai sin(2i//i) a 2 sin(20^/) 

ai ai cos(2j/r/) ai sin(20^/) 0 

, a 2 0^2 cos(2i//i) a 2 sin(20^/) 0 


(14) 


To ameliorate the limited range of orientations, we perform 
a joint solution for all the RCAs at each frequency, in contrast 
to the WMAP approach of solving for each radiometer inde¬ 
pendently. In Fig. 14 we show output maps from the IQUSS 
approach at 30 GHz: Q and U maps (top row); and and 82 
maps (bottom row). Note that Q and U maps are noisier than 
for the nominal mapmaking solution. Over most of the sky the 
resulting maps of spurious signals are still noisy and therefore 
we chose a conservative approach to estimating the a-factor for 
each RCA. This is done with a weighted least-squares fit of the 
leakage map L to the spurious signal Sk(Sk = cikL in the absence 
of errors) using only those pixels with \h\ <15°, since at higher 
Galactic latitudes the foregrounds and hence the spurious sig¬ 
nals are weak and mainly contribute noise to the solutions. Our 
code removes pixels where the condition number for the noise 
covariance matrix is less than a given threshold. This now 
has a negligible effect, since thanks to the modification of the 
Planck scanning strategy after Survey 5, the matrix is very 
well-behaved; even with our conservative limit of 8 x 10“^ for the 
condition number, fewer than 200 pixels are excluded at 44 GHz 
and none at 30 and 70 GHz. Our derived values for the ( 2 -factors 
are listed in Table 14. 

Beam-shape mismatch between radiometer arms can also 
lead to polarization leakage when there are strong intensity gra¬ 
dients. We therefore examined the effect of excluding compact 
sources using the WMAP 7-year point source mask. However, 
this made no significant changes, apart from a dramatically in¬ 
creased uncertainty, so our final values do not use such masking. 
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Fig. 14. IQUSS solution maps at 30 GHz. Top left: Stokes Q. Top right: Stokes U. Bottom left: spurious signal from the first RCA, 
S 1 . Bottom right: spurious signal from the second RCA, 82 - Polarization maps are noisier than the usual mapmaking solution, since 
S 1 and S 2 have to be extracted from the same data. 


Table 14. Bandpass mismatch ^-factors from fitting the leakage 
model map to the spurious maps. 


Horn 

a-factor 

70 GHz 

LFI 18. 

-0.0018 + 0.0022 

LFI 19. 

0.0124 + 0.0024 

LFI 20. 

0.0034 + 0.0024 

LFI 21. 

-0.0115 + 0.0024 

LFI 22. 

0.0039 + 0.0024 

LFI 23. 

0.0057 + 0.0024 

44 GHz 

LFI 24. 

0.0033 + 0.0005 

LFI 25. 

0.0004 + 0.0004 

LFI 26. 

0.0014 + 0.0004 

30 GHz 

LFI 27. 

0.0046 + 0.0002 

LFI 28. 

-0.0089 + 0.0002 


We compared these values with an independent derivation 
based on aperture photometry of bright sources in the IQUSS 
maps, including the Tarantula nebula in the LMC, which lies 
in the “deep” region around the Ecliptic pole, which is scanned 
with multiple different polarimeter orientations across a wide 
range of angles, and hence allows a particularly accurate blind 
separation of spurious signal. Other calibrators were bright H ii 
regions at relatively high Ecliptic latitude, since the range of po¬ 


larization orientations observed increases towards the Ecliptic 
poles. Hii regions were chosen because they have minimal in¬ 
trinsic polarization, but we did not force Q and U to zero in the 
analysis. The f 2 -factors derived from the calibrators were consis¬ 
tent with our preferred values derived from the large area fit, but 
somewhat less precise. 

11.4. Production of correction maps 

The polarization data from a given radiometer constrains one 
Stokes parameter (say Qh) in a frame of reference tied to the 
specific feed horn (or RCA). This is projected onto the sky ac¬ 
cording to the sky orientation of the horn frame. Hence the con¬ 
tribution of the spurious signal from each radiometer is modu¬ 
lated into the Q and U sky pixels by geometric projection factors. 
This modulation can be derived by re-scanning, in a mapmaking 
fashion, the estimated spurious map S = OkL. Instead of an ac¬ 
tual re-scanning, which is time consuming, we create projecting 
maps Aqyu] by solving the mapmaking system 

^Qmip) = (15) 

where are the maps obtained by binning a stream of -1 for 
the “side” and of -rl for the “main” arms separately. Einally the 
correction maps are 

AQ[U] = L X ^ cikAk^QYu^. (16) 

k 

One of the main drawbacks of deriving our L maps at 1° 
resolution emerges at this stage. The correction must be applied 
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to Q and U maps matched in resolution, and so the raw Q and 
U maps are smoothed to give the required 1° FWHM Gaussian 
beam. We can regard the raw maps as the true Q and U sky, 
smoothed with the instrumental beam, plus the leakage term, 
plus noise. The leakage term in the raw maps can be thought of 
as an infinite-resolution leakage sky convolved with the instru¬ 
mental beam, and then multiplied by the leakage projection maps 
Pq[U] = which are defined at the pixel level ac¬ 

cording to Eq. (14). When we smooth this raw map, we smooth 
the product PL, but when we construct our correction map, we 
have only the smoothed L map. The smoothed product is not 
equal to the smoothed L map multiplied by the full-resolution 
P map, which contains fine-scale structure induced by caustics, 
lost data, and abrupt changes in the survey strategy; nor is it 
equal to the smoothed L-maps multiplied by the smoothed P 
map, which is over-smoothed in regions where both L and P 
vary rapidly. In practice we used the smoothed P maps, since 
P only varies rapidly near a small subset of pixels, and only a 
few of these will also have rapidly varying L. The issue is most 
significant for compact sources, for which we recommend anal¬ 
ysis of the raw maps, followed by a leakage correction using the 
derived IQU fluxes, as described in Planck Collaboration XXVI 
(2016). 

Although our ( 2 -factor estimates are relatively stable, our fit 
of the leakage maps to the spurious maps showed significant 
residuals at the level of 18yuK, 24//K, and 16yuK, respectively 
at 30, 44, and 70 GHz. Contributing factors may include: errors 
in the leakage maps caused by errors in the component sepa¬ 
ration; residual beam ellipticity after smoothing to 1°; and any 
variation with time of the bandpass, which would cause corre¬ 
sponding changes in the ^-factors. As a check on our results, 
we used the IQUSP procedure in which we create a prior for 
each S map using our component-separation L map and our best 
( 2 -factor estimates. This process returns the prior 5-map essen¬ 
tially unchanged over most of the sky, and hence gives Q and U 
maps indistinguishable from our corrected versions. However, 
it prefers the IQUSS solution when it differs significantly from 
the prior (essentially in regions of the brightest foreground emis¬ 
sion, where the limitations of our simplified emission models 
become apparent). An advantage of the method is that the maps 
are returned at full resolution, wherever the data can constrain 
the resolution to be higher than in the prior. These maps confirm 
that most details of the structure along the Galactic plane in our 
corrected LFI polarization maps are consistent with the data, i.e. 
are reproduced in the IQUSP images, including the most sig¬ 
nificant discrepancies with WMAP. Although not fully validated 
and therefore not included in the current release, the IQUSP 
images are likely to form the baseline for our final-release polar¬ 
ization maps. 


12. Data validation 

We verify the quality of the LFI data with a suite of null tests, 
as well as with a set of simulations reproducing the main in¬ 
strumental systematic effects and the calibration process. In this 
section we summarize the main results of our analysis, and refer 
to Planck Collaboration III (2016) and Planck Collaboration V 
(2016) for more details. 

Null tests are performed on blocks of data covering differ¬ 
ent time scales (from the pointing period to surveys and years) 
and considering different instrument combinations (radiometer, 
horn, horn-pairs, and frequency) both in total intensity and po¬ 
larization (when applicable). 


Such null tests can probe different systematic effects depend¬ 
ing on the time and instrument selection considered. Differences 
at horn level between odd and even surveys may show effects 
due to the sidelobe contribution, since the relative orientation of 
the horns with respect to the sky is changed. Furthermore, the 
comparison of power spectra at the frequency level may reveal 
the impact of calibration uncertainties related to the relative ori¬ 
entation of the scans and the CMB dipole, our main calibration 
source as discussed in Sect. 7. 


12.1. Null test results 

In order to assess null test results, it is fundamental to define a 
clear figure of merit as a pass-fail criterion. Failure of a spe¬ 
cific test is an indication of a data problem and/or issues in 
data processing that should be studied further. As we already 
did for the previous release, we take the noise level as derived 
from “half-ring” difference maps, made of the first and second 
half of each stable pointing period (half-ring maps) weighted by 
the hit count, as the figure of merit. This quantity traces the ac¬ 
tual properties of the data, including white noise, as well as un¬ 
modelled and un-corrected effects. Figure 15 shows results at 
the frequency level for both TT and EE power spectra when we 
compare survey differences to the noise level derived from the 
corresponding half-ring maps. For simplicity we show here only 
a subset of survey differences that are illustrative of the general 
trend. 

When interpreting these results it is important to note that 
we have substantially improved the quality of data at 30 GHz 
by using the new An calibration (Planck Collaboration V 2016), 
which accounts for the impact of the full beam during calibra¬ 
tion, and by removing at the TOD level the modelled sidelobe 
signals of both the CMB dipole and Galactic emission (derived 
from the FFP8 simulation runs). This is particularly evident from 
TT spectra, where the null test data match the level of the half¬ 
ring differences. However, there is still an issue in polariza¬ 
tion when considering differences involving Surveys 2 and 4. 
At 44 GHz, which has the lowest sidelobes among LFI chan¬ 
nels, the agreement with the half-ring noise is extremely good. 
We have almost the same situation at 70 GHz, although at very 
low multipoles {i < 10) there are discrepancies between survey 
difference and half-ring noise; again this is particularly evident 
when Surveys 2 and 4 are considered. 

To be more quantitative about these results at 70 GHz we 
compute deviations from the half-ring noise in terms of 




(17) 


We sum up single values in the range 2-50 and from;^2 and 
Adof we derive p-values of the distribution. While in principle 
proper noise simulations should be used, for our purposes it is 
sufficient to consider simple half-ring noise, which is already 
able to reveal interesting features in the data. Table 15 reports 
and p-values for the three survey differences shown in Fig. 15, 
which suggests that Surveys 2 and 4 clearly yield poor and 
problematic p-values. 

As discussed further in Sect. 13, on the basis of these and 
other results, we have discarded these two surveys from the re¬ 
leased likelihood. 
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Fig. 15. Null test results comparing power spectra from survey differences to those from the half-ring maps. Differences are: left 
Survey 1 - Survey 2; middle Survey 1 - Survey 3; and right Survey 1 - Survey 4. These are shown for 30 GHz (top), 44 GHz 
(middle), and 70 GHz (bottom) for both TT and EE power spectra. 


Table 15. Survey differenceand p-values. 


Survey differences 


p-value 

Survey 1 - Survey 2 

83.3 

0.0016 

Survey 1 - Survey 3 

70.8 

0.0221 

Survey 1 - Survey 4 

98.3 

3.7 X 10-5 


12.2. Half-ring test 

As already pointed out, the half-ring difference maps are the 
best direct information about the actual noise in the LFI data. A 
proper characterization of the noise is fundamental for the cre¬ 
ation of realistic NCVMs and noise MC required for the follow¬ 
ing steps in the data analysis. In this respect such noise mod¬ 
elling has to be validated against the half-ring maps. For the cur¬ 
rent analysis we followed the same procedure exploited in the 
previous data release. We computed auto-spectra in temperature 
and polarization with anafast of both the half-ring difference 
maps and 10000 noise Monte Carlo simulated maps taken from 
FFP8. We compared the half-ring spectra with the distribution 
of the noise MC simulations and with the white noise derived 


from the white noise covariance matrices (WNCVM) calculated 
by Madam map-making. 

Figure 16 gives a flavour of this comparison for the three 
LFI frequencies and for both total intensity TT and polarization 
(EE and BE) power spectra. Note that the half-ring noise spectra 
are binned over a range of A^ = 25 for ^ > 75. The agreement 
between half-ring noise spectra and noise MC distribution is re¬ 
markable, and gives us confidence that the LFI noise properties 
are accurately characterized. 

We further inspect this comparison computing the mean Q 
for the high-^ tail of the spectrum (1150 < ^ < 1800) and com¬ 
paring it with the WNCVM (white noise covariance matrix) es¬ 
timate (Fig. 17). It is clear that there is some residual 1// noise 
also at high-^ as has been already pointed out in the 2013 re¬ 
lease (Planck Collaboration II 2014). This means that both data 
and noise MCs predict a slightly higher noise than the WNCVM. 
The residual is of the order of 1.6 % (TT) maximum at 30 GHz, 
1.3 % (BB) at 44 GHz, and 1.0 % (EE) at 70 GHz. On the other 
hand, the agreement between the actual data and the full noise 
MCs is extremely good being of the order of 0.5 % at 30 GHz, 
0.4 % at 44 GHz, and 0.2 % at 70 GHz. 
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Fig. 16. Consistency of the noise angular power spectra from the half-ring difference maps (red), white noise covariance matrix 
(black dash-dotted line), and 10000 full-noise MCs (grey band showing 50% quantiles, black solid line, and limits at 16% and 
84 % quantiles). From top to bottom we have TT, EE, and BB spectra for 30 (left), 44 (centre) and 70 GHz (right). Half-ring spectra 
are binned with = 75 for ^ > 75. 



Frequency channel 


Fig. 17. Ratio of the mean noise angular power spectrum in the 
high-T” (1150 < ^ < 1800) tail to the white noise as derived from 
the white noise covariance matrices from Madam. 


12.3. End-to-end test results 

The LFI calibration pipeline is necessarily quite complex, since 
it includes iterative mapmaking, sidelobe removal. Galactic 
masking, map domain fits to the 471 beam-convolved dipole, and 
filtering. While the accuracy of the mean calibration constant is 
important, particularly for inter-frequency validation and fore¬ 
ground modelling, we are mostly concerned with quantifying the 
level of systematic errors in our estimation of the calibration over 
time. The gain of the LFI radiometers typically varied by a few 
percent over the four years mission lifetime, with changes at time 
scales from single pointing periods to the full mission. Null tests 
on survey and year time scales set useful limits on systematic 
effects, including incorrect calibration estimation, but it is still 
important to develop a “bottom up” estimation of possible er¬ 
rors. Consequently we have carried out several parallel efforts to 
simulate our calibration procedure, each using different software 
and detailed choices for inputs, but following the same general 
approach, which we now summarize. 


- Start with a fiducial sky map (in kelvins), either from the 
frequency maps of the data, FFP8, or some other simula¬ 
tion. This map includes CMB anisotropies, foregrounds, and 
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possibly some systematics, but no dipole signals, and can be 
either temperature only or Q and U. 

- “Unwrap” or rescan the map to a time ordered signal data 
set, in “ring” basis (still in Kelvin). This is done using actual 
flight pointing data. 

- Add dipole signals, including the solar dipole and the orbital 
dipole. We can choose here whether to use a “pencil beam” 
model, where the dipole signal that is added has been sam¬ 
pled from the sky model with a Dirac delta function, or a 
An model consisting of an all-sky convolution of the detector 
beam model with the dipole model. 

- Add instrument noise, either white noise or full 1// noise 
(only white noise turns out to be relevant). 

- All the steps described so far assume a timeline in kelvins. 
Next we “decalibrate” these simulated data streams using a 
fiducial model for the actual detector gain, and produce time¬ 
lines in volts. A standard choice here is to use the so-called 
“Delta V” gain, which is a radiometer gain estimated directly 
from the DC-coupled detector data. While we know from de¬ 
tailed tests that this gain does not track the actual gain fiuctu- 
ations better than about 0.5 %, it has the advantage of being a 
gain estimate with no smoothing applied, and should reflect 
closely the true statistics of the radiometer gain. 

- From this simulated timestream, we proceed with our nomi¬ 
nal calibration pipeline to recover the input gain. In this way, 
we can compare the recovered time domain gain estimate 
to the fiducial input, as well as the final calibrated maps to 
the fiducial input maps. The results of such comparisons are 
shown for two radiometers in Figs. 18 and 19. 


These simulations are designed to test the impact of our pro¬ 
cedures on the results. They also provide a mechanism for quan¬ 
titatively determining the impact that errors in the inputs, such 
as beam shape or far sidelobe contribution, have on our output 
maps and other scientific products. They do not provide a way to 
estimate what those input errors are; these must be determined 
by dedicated investigations on the optical model or instrument- 
specific simulations. Starting from reasonable estimates of the 
systematics affecting our instrument, however, we can introduce 
changes in the input within the expected range and then test for 
deviations in the recovered calibration. We thus obtain both the 
sensitivity to that effect and an estimate of the probable error 
causing it, assuming either extreme values (conservative) or the 
expected 1 cr (typical). Similarly, we can use this approach to 
determine the sensitivity of the calibration process to Galactic 
masking. 

The basic results of such end-to-end tests of the effects of 
systematics are summarized in Table 16. Comparison of the dif¬ 
ference between input and output gains shows a typical bias of 
order 0.2 %. 

Figures 18 and 19 show input and output gain constants and 
the relative variations between the two, at 30 GHz and 70 GHz, 
respectively. The 30 GHz channels are the most difficult to cal¬ 
ibrate, because they are more sensitive to changes in instru¬ 
ment configuration, causing a bigger number of jumps, while 
the 70 GHz are more sensitive to the instrument noise. Table 16 
shows mean and standard deviations of the relative variations 
between the input gain constants and the output ones, for all the 
radiometers. The resulting precision of the photometric calibra¬ 
tion is up to 0.2 %, thus validating the calibration algorithm. 


12.4. Intra-frequency consistency check 

We tested consistency between 30, 44, and 70 GHz maps by 
means of power spectra, as already done in the previous release 
(Planck Collaboration III 2014). In order to avoid the need to es¬ 
timate the noise bias, we simply took the cross-spectra between 
half-ring maps at the three LFI frequencies. As in the 2013 data 
release, we used the cROMAster code which extends the pseudo- 
Q approach of Hivon et al. (2002) to cross-power spectrum es¬ 
timation (Polenta et al. 2005). Although suboptimal with respect 
to the maximum likelihood approach, this method provides ac¬ 
curate results, and is at the same time computationally quick 
and light. Consequently, this method is widely used within the 
CMB community (see e.g. Molinari et al. (2014) and references 
therein for a comparison between different power spectrum esti¬ 
mators). 



Pid 


Fig. 18. Relative variations between input and output of the end- 
to-end test for radiometer 27S at 30 GHz. In general, we recover- 
the input gain to better than 0.1 %, except for some larger excur¬ 
sions introduced by sudden changes in the instrument configura¬ 
tion, to which the 30 GHz radiometers are particularly sensitive. 



Pid 


Fig. 19. Relative variations between input and output of the end- 
to-end test for radiometer 22S at 70 GHz. The overall recovery 
is under 0.1 %, with some spikes in the longest pointing periods. 
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Table 16. Mean and associated error of the percentage variation 
between input and output of the end to end tests 


Radiometer 

Mean difference [%] 

70 GHz 

LFI 18M 

-0.002 + 0.057 

LFI 18S 

0.040 + 0.140 

LFI 19M 

0.007 + 0.169 

LFI 19S 

0.058 + 0.090 

LFI 20M 

0.009 + 0.081 

LFI 20S 

-0.002 + 0.071 

LFI 21M 

-0.005 + 0.070 

LFI21S 

-0.031 + 0.066 

LFI 22M 

0.012 + 0.093 

LFI 22S 

0.007 + 0.062 

LFI 23M 

0.016 + 0.087 

LFI23S 

0.029 + 0.083 

44 GHz 

LFI 24M 

-0.001 + 0.098 

LFI 24S 

-0.023 + 0.079 

LFI 25M 

-0.009 + 0.074 

LFI 25S 

0.015 + 0.098 

LFI 26M 

0.034 + 0.084 

LFI 26S 

-0.027 + 0.082 

30 GHz 

LFI 27M 

0.040 + 0.094 

LFI 27S 

-0.001 + 0.127 

LFI 28M 

-0.034 + 0.092 

LFI 28S 

0.052 + 0.136 


Those spectra are computed using a mask that is the com¬ 
bination of the G040, G060 and GOTO Planck masks, respec¬ 
tively at 30, 44, and 70 GHz, together with the proper frequency- 
dependent point source mask. In Fig. 20 cross-spectra from 30, 
44, and 70 GHz half-ring maps are presented, showing very good 
agreement among these maps (especially as we did not apply 
any component separation to the maps). All three data sets show 
strong consistency with the Planck best-fit TT spectrum (black 
points) to which a contribution from unmasked point sources has 
been properly added. 

Another more quantitative way for assessing data consis¬ 
tency is to build scatter, or TT-, plots for the three frequency 
pairs. In order to do this we have to subtract the contribution 
of point sources below the mask threshold at each individual 
frequency. After that we perform a linear fit, accounting for er¬ 
rors in both X- and y-axes, to quantify the level of agreement 
between pairs. Results are presented in Fig. 21, where we com¬ 
pare spectra in the multipole range around the first acoustic peak. 
The agreement is extremely good and spectra are consistent with 
unity within the errors (deviations are between 0.9 and 0.1 %). 
That in turn means a calibration accuracy in the map at the sub¬ 
percent level. This is very significant considering that we did 
not take into account foreground removal or uncertainties on the 
window function and calibration; therefore we may expect the 
agreement to improve when these issues are taken into account. 



Fig. 20. Temperature cross-power spectra (from half-ring maps) 
at 30, 44, and 70 GHz, binned in multipole space. Foreground 
emission is excluded only by means of a Galactic sky mask, 
without further component separation. Best-fit Planck temper¬ 
ature spectra plus contributions from un-resolved point sources 
are shown as dashed lines for each LFI band. 


12.5. Internal consistency check 

In order to assess the internal consistency of 70 GHz data, we 
build three flavours of cross-power spectra that use different 
kind of data splits, namely the half-ring maps, the detector set 
(quadruplet) maps, and the year 1-3 and year 2-4 maps. In 
Fig. 22 we show residuals of the three estimates compared to 
the expected deviations computed by running the same proce¬ 
dure on the realistic FFP8 Monte Carlo simulations. A simple 

analysis shows that residuals are compatible with the null hy¬ 
pothesis. 

We then apply the Hausman test (Polenta et al. 2005) to fur¬ 
ther verify the consistency of the three power-spectrum esti¬ 
mates. We define the statistic 


Ht = {C( - Ce) / ^Var{Q-C(}, 


(18) 


where Q and Q represent two different cross-spectra, and we 
combine the information from different multipoles through the 
quantity 



(19) 


where [.] denotes integer part. It can be shown that the distribu¬ 
tion of Biir) converges to a Brownian motion process, which can 
be studied using three test statistics defined as = sup 

S 2 = sup^|5z,(r)| and S 2 , = B^(r)dr. Results for the compari¬ 
son of detector set (DS) and year based (YR) cross-spectra are 
shown in Fig. 23. Vertical lines represent the values of the test 
statistics computed from Planck maps as compared to the em¬ 
pirical distribution of the test statistics derived from FFP8 sim¬ 
ulations. The application of the Hausman test to the other cross¬ 
spectra combinations produces similar results, thus supporting 
the strong internal consistency of the LFI 70 GHz data. 

In this second Planck data release, the calibration pipeline 
considers the full convolution between the beam response B and 
the calibration signal D. This is a novel approach, which allows 
us to better control the impact of optical systematic effects on 
the calibration and to improves the self-consistency of the data. 
Note that in the first data release, the dipole fitting routines used 
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Fig. 21. Consistency between cross-power spectra at LFI frequencies: left 70 GHz versus 30 GHz; middle 70 GHz versus 44 GHz; 
and right 44 GHz versus 30 GHz. The solid red line is the linear regression, accounting for error on both axes. Slope values are 
found to be consistent within the uncertainties. 



Fig. 22. Residuals between three different cross-power spectra 
computed from 70 GHz data: half-ring (HR) maps, quadruplet 
(detector set, DS) maps, and year 1-3 / year 2-4 (YR) maps. Error 
bars are derived from the realistic FFP8 simulations. 


to measure the calibration constants assumed a pencil-like beam, 
and the mismatch in power was fixed by applying a beam win¬ 
dow function to the power spectra. 

As Planck Collaboration V (2014) has shown, the convolu¬ 
tion B ^ D retains the same dipole shape as D, but there are two 
effects of particular relevance for this discussion: 

1. the finite width of the main beam and the presence of lobes 
reduces the peak-to-peak amplitude of the dipole itself (i.e. 
the peak-to-peak variation in 5 * D is smaller than the varia¬ 
tion in D); 

2. the lack of perfect axial symmetry (particularly in the region 
which is far from the main beam) induces a tilt in the dipole 
axis. 


The first point implies that using the B^D signal as a calibration 
source reduces the average value of the calibration constant K 
([K] = KV-i). Planck Collaboration V (2014, 2016) quantify 
the amount of such variation in terms of the measured power 
spectra Q: 

^ 1 ~ /si - 0sky + 0D 1 

f(pencn)-[ J’ 

where /si is the fraction of B that falls outside 5° of the main 
beam (the “sidelobes”), fo = dtB^\ * T>/5?5main * T) is the ratio 
between the variation of the dipole signal entering the sidelobes 
and the variation of the same signal entering the main beam, 
and /sky and are defined similarly io fo but in terms of the 
amount of Galactic signal plus CMB (/sky), and of the CMB 
alone 

We have verified the consistency of this approach by pro¬ 
ducing a set of maps using data from the current release, but 
calibrated using the pencil-beam approximation. By comparing 
the raw power spectra of these maps with the official LFI power 
spectra of the second release, we have measured excellent agree¬ 
ment (better than 0.03 %) with the estimate provided by Eq. (20), 
apart from four out of six 44 GHz radiometers. In the 44 GHz 
case, however, because of the small level of the sidelobes, the 
resulting change in the Q is (at < 0.4 %) still smaller than for 
the other two LFI bands. 

12.6. Updated systematic effects assessment 

Known instrumental systematics affecting LFI maps are dis¬ 
cussed in detail in Planck Collaboration III (2016) and are listed 
in Table 17, along with short descriptions of their causes and 
strategies for their removal. In Tables 18, 19, and 20 we list both 
the rms and the difference between the 99 % and the 1 % quan¬ 
tiles in the pixel value distribution for the /, Q, and U maps, 
at 30, 44 and 70 GHz respectively. We refer to the latter as the 
peak-to-peak (p-p) difference, even though it neglects outliers, 
since it effectively approximates the peak-to-peak variation of 
the effect on the map. 

Detailed analysis reported in Planck Collaboration III (2016) 
shows that systematic uncertainties are at least two orders of 
magnitude below the CMB TT power spectrum and are not sig¬ 
nificantly contaminating the EE and BB spectra. 


^ Refer to Planck Collaboration V (2016) for a mathematical deriva¬ 
tion and a discussion of the formula. 
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Fig. 23. From Left to Right, the empirical distribution (estimated via FFP 8 simulations) of the ^' 2 , ^3 statistics of the Hausman 
test (see text). Vertical lines represent the values obtained from Planck 70 GHz data. 

Table 17. List of known instrumental systematic effects in Planck-h¥l. 


Effect 

Source 

Control/removal 


Effects independent of sky signal (T and P) 

White noise correlation 

Phase switch imbalance 

Diode weighting 

1// noise 

RF amplifiers 

Pseudo-correlation and destriping 

Bias fluctuations 

RF amplifiers, back-end electronics 

Pseudo-correlation and destriping 

Thermal fluctuations 

4-K, 20-K and 300-K thermal stages 

Calibration, destriping 

1-Hz spikes 

Back-end electronics 

Template fitting and removal 


Effects dependent on the sky signal (T and P) 

Main beam ellipticity 

Main beams 

Accounted for in window function 

Near sidelobe 
pickup 

Optical response at angles 
<5° from the main beam 

Masking of Galaxy and point 
sources 

Far sidelobe pickup 

Main and sub-reflector spillover 

Model sidelobes removed from timelines 

Analogue-to-digital 
converter non-linearity 

Back-end analogue-to-digital 
converter 

Template fitting and removal 

Imperfect photometric 
calibration 

Sidelobe pickup, radiometer noise 
temperature changes and other 
non-idealities 

Adaptive smoothing algorithm using An 
beam, 4-K reference load voltage output 
temperature sensor data 

Pointing 

Uncertainties in pointing reconstru¬ 
ction, thermal changes affecting 
focal plane geometry 

Negligible impact on anisotropy 
measurements 


Effects specifically impacting polarization 

Bandpass asymmetries 

Differential orthomode transducer 
and receiver bandpass response 

Spurious polarization removal 

Polarization angle 
uncertainty 

Uncertainty in the polarization 
angle in-flight measurement 

Negligible impact 

Orthomode transducer 
cross-polarization 

Imperfect polarization separation 

Negligible impact 


13. Low-^ data selection 

The 70 GHz polarization data are of special importance since the 
Planck low-^ likelihood (Planck Collaboration XI 2016) used to 
determine cosmological parameters is based on them. In order to 
provide the best data possible for the construction of the likeli¬ 
hood, we perform several tests at survey level in order to choose 
the most reliable data combination. For this purpose we focus on 
the very low multipoles, especially £ = 2-4, which are the most 
susceptible to systematic errors. 


We compare results from actual data and from noise- 
only Monte Carlo realizations made for the FFP 8 simulations. 
Specifically, we take differences between the full data set (over 
the entire mission lifetime) and some specific combinations of 
surveys, for both noise simulations and real data. We then com¬ 
pute the angular power spectra of these differences to look for 
anomalies. 

The analysis at the level of surveys is very informative: as 
a consequence of the scanning strategy and payload geome¬ 
try, Survey 1 and Survey 3 share the same beam orientation 
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Table 18. Summary of systematic effect uncertainties on 30 GHz 
maps^ in //Kcmb- Columns give the peak-to-peak (“p-p”) and 
rms levels for Stokes /, 2, and U maps. 



1 

r 


Q 


U 


p-p 

rms 

p-p 

rms 

p-p 

rms 

Near sidelobes. 

. 0.72 

0.13 

0.05 

0.01 

0.05 

0.01 

Pointing . 

. 0.37 

0.07 

0.02 

0.01 

0.02 

0.00 

Polarization angle . . . 

. 0.02 

0.00 

0.53 

0.11 

0.64 

0.15 

1-Hz spikes. 

. 0.54 

0.11 

0.11 

0.02 

0.09 

0.02 

Bias fluctuations .... 

. 0.07 

0.01 

0.07 

0.01 

0.06 

0.01 

ADC non-linearity . . . 

. 0.42 

0.09 

0.54 

0.11 

0.56 

0.11 

Calibration. 

. 2.43 

0.55 

2.53 

0.46 

2.34 

0.43 

Thermal fluct. (300 K) 

. 0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Thermal fluct. (20 K) . 

. 0.12 

0.03 

0.06 

0.02 

0.06 

0.02 

Thermal fluct. (4K) . . 

. 0.29 

0.06 

0.06 

0.01 

0.05 

0.01 

TotaP. 

. 2.72 

0.61 

2.79 

0.52 

2.42 

0.49 


^ Calculated for a pixel size approximately equal to the average beam 
FWHM. A null value indicates a residual < 10"^ ^uKcmb- 
^ The total has been computed on maps resulting from the sum of indi¬ 
vidual systematic effect maps. 

Table 19. Summary of systematic effect uncertainties on 44 GHz 
maps in /iKcMB- Columns give the peak-to-peak (“p-p”) and rms 
levels for Stokes /, 2, and U maps. 




I 

Q 


U 


p-p 

rms 

P-P 

rms 

p-p 

rms 

Near sidelobes. 

. 0.09 

0.02 

0.00 

0.00 

0.00 

0.00 

Pointing . 

. 0.30 

0.06 

0.01 

0.00 

0.01 

0.00 

Polarization angle . . . 

. 0.04 

0.01 

0.35 

0.07 

0.38 

0.10 

1-Hz spikes. 

. 1.99 

0.40 

0.88 

0.18 

1.04 

0.21 

Bias fluctuations .... 

. 0.04 

0.01 

0.05 

0.01 

0.05 

0.01 

ADC non-linearity . . . 

. 0.30 

0.06 

0.36 

0.07 

0.34 

0.07 

Calibration. 

. 1.05 

0.18 

1.57 

0.29 

1.31 

0.26 

Thermal fluct. (300 K) 

. 0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Thermal fluct. (20 K) . 

. 0.04 

0.02 

0.06 

0.01 

0.05 

0.01 

Thermal fluct. (4 K) . . 

. 0.23 

0.05 

0.05 

0.01 

0.06 

0.01 

Total. 

. 2.29 

0.45 

1.95 

0.37 

1.76 

0.37 


with respect to the sky. The same is true for Surveys 2/4, 5/7, 
and 6/8. For this reason we consider these combinations jointly 
for the null tests, thus maximizing signal-to-noise. Figure 24 
shows the distribution of angular power for E- and 5-modes for 
each survey pair, as derived from the Monte Carlo simulations, 
with results from the actual LFI data indicated by vertical lines. 
Evidently, Survey 2 and Survey 4 are quite anomalous with re¬ 
spect to the rest of the surveys. We will offer some possible ex¬ 
planations below. First, however, we can be more quantitative 
and compute the probability to exceed (PTE) of our data, based 


Table 20. Summary of systematic effect uncertainties on 70 GHz 
maps in //Kcmb- Columns give the peak-to-peak (“p-p”) and rms 
levels for Stokes /, Q, and U maps. 



i 

r 


Q 


U 


p-p 

rms 

p-p 

rms 

p-p 

rms 

Near sidelobes. 

. 0.30 

0.07 

0.01 

0.00 

0.01 

0.00 

Pointing . 

. 0.60 

0.11 

0.03 

0.01 

0.03 

0.01 

Polarization angle . . . 

. 0.02 

0.00 

0.08 

0.02 

0.08 

0.02 

1-Hz spikes. 

. 0.39 

0.08 

0.17 

0.03 

0.15 

0.03 

Bias fluctuations .... 

. 0.68 

0.14 

0.84 

0.17 

0.95 

0.18 

ADC non-linearity . . . 

. 1.56 

0.33 

1.92 

0.39 

2.05 

0.41 

Calibration. 

. 1.06 

0.23 

0.98 

0.18 

0.77 

0.16 

Thermal fluct. (300 K) 

. 0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Thermal fluct. (20 K) . 

. 0.44 

0.08 

0.07 

0.01 

0.08 

0.02 

Thermal fluct. (4 K) . . 

. 0.38 

0.08 

0.04 

0.01 

0.05 

0.01 

Total. 

. 2.24 

0.47 

2.27 

0.46 

2.38 

0.48 


on simulations. Results are reported in Table 21. These probabil¬ 
ity values seem to indicate that Surveys 2 and 4 show systematic 
effects. Guided by these findings, we report the PTE values for 
the differences between the full mission and the survey combi¬ 
nations in Table 22. 

We can also combine the PTE results from the survey null 
tests across these multipoles. In Table 23 we report results from 
a test of uniformity of the PTEs, simply counting how many en¬ 
tries are lower than a given threshold. The -values for these 
tests are computed assuming binomial statistics. We report re¬ 
sults for different values of the threshold to show their robustness 
and stability with respect to the thresholds. 

Quantitatively Surveys 2 and 4 again stand out as anoma¬ 
lous at roughly the 3 a level. Currently the reason for this is not 
fully understood, but we note that this particular survey pair has 
a scanning strategy that produces larger uncertainties in gain, as 
demonstrated in Fig. 4. The geometry for these two surveys also 
increases the sensitivity of the very low-^ results to small er¬ 
rors in estimates of Galactic contamination of the far sidelobes. 
These issues are under investigation and will be addressed fur¬ 
ther in the next data release, but for the moment we choose to be 
conservative and remove Surveys 2 and 4 from the default likeli¬ 
hood developed in Planck Collaboration XI (2016). The default 
likelihood is used in Planck Collaboration XIII (2016) to derive 
cosmological parameters. The optical depth to reionization, r, is 
the parameter most affected by this choice: removing Surveys 2 
and 4 changes the value of this parameter by about 0.5 cr. 

14. The low-7^ likelihood 

The baseline 2015 Planck low-^ likelihood is described in depth 
in Planck Collaboration XI (2016). Here we briefly discuss its 
polarization content, largely based on data from the Planck 
70 GHz channel. As noted in the previous sections. Survey 2 
and 4 are excerpted from the data set to reduce the chance 
of systematic contamination. In this section, we do not focus 
on the low-^ temperature block of the likelihood developed in 
Planck Collaboration XI (2016) that is based on a CMB map de- 
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Fig. 24. Measured LFI 70 GHz EE {top) and BB {bottom) null power spectra for 7’ = 2, 3, and 4 (vertical lines), compared to 
the distribution derived from noise-only Monte Carlo simulations. Null spectra from the difference between full data and specific 
surveys combinations: left Survey 1 and Survey 3; {middle) Survey 2 and Survey 4; and {right) Survey 5 and Survey 7. It is clear 
that Survey 2/Survey 4 stands out with respect to the others. 


Table 21. PTE for EE and BB low multipoles, for the differences between full mission and individual surveys. 


Multipole 

Full - S 1 

Full -S 3 

Full - S 2 

Full - 5 4 

Full - S 5 

Full — S'] 

Full - 5 6 

Full - 5 8 

EE 

. 

0.885 

0.307 

0.328 

0.015 

0.241 

0.975 

0.837 

0.090 

£ = 3 . 

0.740 

0.730 

0.137 

0.223 

0.206 

0.566 

0.377 

0.064 

£ = 4 . 

0.807 

0.828 

0.890 

0.535 

0.290 

0.998 

0.932 

0.476 

BB 

. 

0.998 

0.214 

0.030 

0.482 

0.098 

0.680 

0.986 

0.092 

^ = 3. 

0.833 

0.796 

0.843 

0.002 

0.414 

0.823 

0.516 

0.255 

^ = 4. 

0.903 

0.105 

0.092 

0.399 

0.052 

0.524 

0.862 

0.950 


Table 22. PTE for EE and BB low multipoles, for the differences between full mission and survey combinations. 



Multipole 

Full - {Si + ^ 3 ) 

Full — {S 2 + 5 * 4 ) 

Full — {S 3 + S ']) 

Full — {S 6 “1“ 5* g) 

EE 


£^2 . 

0.491 

0.114 

0.526 

0.578 


. 

0.528 

0.137 

0.109 

0.598 


^ = 4. 

0.750 

0.825 

0.976 

0.978 

BB 


£^2 . 

0.482 

0.023 

0.156 

0.544 


. 

0.698 

0.010 

0.866 

0.320 


^ = 4. 

0.218 

0.190 

0.152 

0.995 


rived using the Commander algorithm which employs all Planck At multipoles i < 30, we model the likelihood assuming 
channels from 30 to 353 GHz (Planck Collaboration IX 2016). that the maps are Gaussian distributed with known covariance 

(Planck Collaboration XV 2014): 


£{€,) = P{m\C,) = 


1 

27rV2|M|i/2 



( 21 ) 


31 































































































Planck Collaboration: LFI data processing 


Table 23. Uniformity of the PTEs for survey null tests based on 
the number of entries lower than a given threshold (/^-values are 
from the binomial distribution). 


Threshold 

SS1/SS3 

SS2/SS4 

SS5/SS7 

SS6/SS8 

A < 0.02 

0 (0.215) 

2 (0.002) 

0 (0.215) 

0 (0.215) 

A < 0.05 

0 (0.456) 

3 (0.002) 

0 (0.456) 

0 (0.456) 

A < 0.10 

0 (0.716) 

4 (0.004) 

2 (0.111) 

3 (0.026) 

A < 0.25 

2 (0.609) 

6 (0.014) 

4(0.158) 

3 (0.351) 


where n is the total number of observed pixels, M(Q) is the co- 
variance matrix of m = [T, Q, U], being T, Q, and U the pixel 
space intensity and linear polarization Stokes parameter maps. 
Note that the covariance matrix depends on the CMB model an¬ 
gular power spectra, Q, only through the CMB signal covari¬ 
ance matrix: 

M(C,) = S(C^) + N . (22) 

In order to clean the 70 GHz Q and U maps, we per¬ 
form a template fitting procedure using the Planck 30 GHz 
channel as a tracer of polarized synchrotron emission and the 
Planck 353 GHz channel as a tracer of polarized dust emission. 
Restricting from now onwards m to the Q and U maps (i.e. 
m = [ 2 , f/]) we write: 

m = -—- (mvo - am^Q - , (23) 

1 - a- p 

where mvo, m^o and are bandpass corrected versions 

of the 70, 30, 353 maps (Planck Collaboration III 2016; 
Planck Collaboration VII 2016), whereas a and p are the scal¬ 
ing coefficients for synchrotron and dust emission, respectively. 
The latter are best fitted by minimizing the quantity 

= {\-a-Pfm^ [S(Q) + Nvo]'* m, (24) 

where N 70 is the pure polarization part of the 70 GHz 
noise covariance matrix^ (Planck Collaboration VI 2016), 
and Cp is taken as the Planck 2013 fiducial model 
(Planck Collaboration XVI 2014). We have verified that chang¬ 
ing this model does not impact the results significantly. We 
find a = 0.063, p = 0.0077, with three sigma uncertainties 
(Ta = 0.025 and Oyj = 0.0022. The best fit values quoted corre¬ 
spond to a polarization mask that allows 47% of the sky to pass 
through. In fact, we have repeated this procedure for a set of 24 
masks, allowing sky fractions from 80% to 29%. Such masks 
have been constructed by rescaling the templates and 
to 70 GHz assuming fiducial spectral indexes, computing the po¬ 
larized intensity P = ^JQp + U^ and thresholding the latter. For 
each mask, we evaluate the probability to exceed > x^)' 
The 47% analysis mask is chosen as the tightest mask satisfying 
the requirement !P > 5%. 

We define the final polarization noise covariance matrix used 
in Eq. (22) as: 

N = _ 0^.2 (^™ + Cr^OT353m353) . (25) 

t i d p) 

We have verified that the external (column to row) products in¬ 
volving the foreground templates are subdominant corrections. 
We do not include further correction terms associated with the 
band pass leakage error budget since they are completely negli¬ 
gible. 

^ We assume that the noise induced TQ and TU correlations are neg¬ 
ligible. 


15. Product description 

We now give a list and brief description of Planck LFI released 
products, which can be freely accessed via the Planck Legacy 
Archive interface^, based on all the data acquired during rou¬ 
tine operation from 12 August 2009 to 23 October 2014; the full 
format is reported in the Explanatory Supplement^. 

- Pointing timelines: one FITS file for each OD for each fre¬ 
quency, each FITS file contains the OBT (onboard time) and 
the three angles, 0 , 0, and 0, which identify each sample on 
the sky. 

- Time timelines: one FITS file for each OD for each fre¬ 
quency, each FITS files containing the OBT and its corre¬ 
sponding TAI (International Atomic Time) value (with no 
leap second) in modified Julian day format. This will allow 
the user to cross-correlate OBT with UTC. 

- Housekeeping timelines: all the housekeeping parameters 
with their raw and calibrated values are provided, separated 
by the housekeeping sources and for each OD. 

- Timelines in volts: raw scientific data in engineering units for 
each detector at 30, 44, and 70 GHz and each OD, before its 
calibration from which instrumental systematic effects have 
been removed; 

- Cleaned and calibrated timelines: provided in Kcmb for each 
detector at 30, 44, and 70 GHz and each OD, after scientific 
calibration from which the convolved dipole and convolved 
Galactic straylight have been removed. 

- Scanning beam: An beam representation used in the calibra¬ 
tion pipeline. 

- Effective beam: sky beam representation as a projection of 
the scanning beam on the maps. 

- Full sky maps at each frequency: maps of the sky at 30, 44, 
and 70 GHz in temperature and polarization at Aside = 1024 
and, in the case of 70 GHz at Aside = 2048. Maps are pro¬ 
vided for different data periods, as detailed in Table 13. Note 
that the polarization convention used for the Planck maps 
is referred to as “COSMO” instead of the “lAU”; see the 
Explanatory Supplement for details. 

- Baseline timelines for the full and half-ring periods: these 
timelines have the baseline offset removed (the length is 
specified in Table 12) during the mapmaking process. 

- Low-resolution maps: maps provided at Aside = 16 and their 
associated full noise covariance matrices. 

- RIMO (reduced instrument model): model provided with all 
parameters that identify the main instrument characteristics 
from noise to bandpass and beam function. 

16. Discussion and conciusions 

We have summarized in this paper all the steps taken to 
assemble, calibrate, and map the data gathered by the 
Planck LFI instrument. While the focus is on the changes 
in data and methods since our previous release in 2013 
(Planck Collaboration II 2014), this paper provides a complete, 
if brief, description of LFI data processing, and of the resulting 
temperature and polarization maps at 30, 44, and 70 GHz. 
Many supporting details are provided in four additional papers 
accompanying this release, Planck Collaboration III (2016), 
Planck Collaboration IV (2016), Planck Collaboration V 
(2016), and Planck Collaboration VI (2016), which treat 

^ http://archive.esac.esa.int/pla2 

^ http://pla.esac.esa.int/pla/index.html 
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systematic effects, beams, calibration and mapmaking, re¬ 
spectively. We note that Planck Collaboration VII (2016) and 
Planck Collaboration VIII (2016) cover the same set of topics 
for the Planck HFI instrument. 

16.1. Operations, TOI, and beams 

LFI operated stably for all four years of observations (eight sky 
surveys). The last four surveys were performed with a different 
phase angle (see Sect. 2), allowing us to investigate some sys¬ 
tematic effects (and also reducing Galactic stray light). The most 
significant change in LFI operations was the gradual degradation 
of the sorption cooler, and its replacement by a second cooler 
(on OD 460). For the current release, construction of the satel¬ 
lite attitude and pointing takes account of two additional vari¬ 
ables, the distance to the Sun and the temperature of the REBA 
(Planck Collaboration I 2016). 

Routine spacecraft manoeuvers made approximately 8 % of 
the data unusable; other losses of TOI data were < 1 % for all 
three LFI bands. 

The TOI required several small corrections described in 
Sect. 4. These include corrections for ADC non-linearity and for 
electronic spikes. Residual effects in the LFI maps are at the pK 
level or below (see Planck Collaboration IV (2016) for a fuller 
discussion). 

Measurements of LFI beam properties (Sect. 6) have sub¬ 
stantially improved since the earlier release, based on repeated 
scans of Jupiter and better modelling of sidelobes. The effective 
beam solid angles at 30, 44, and 70 GHz are 1190.06, 832.00, 
and 200.90 [arcmin^], respectively see Tab. 7 for details. The 
remaining sidelobe power outside the main beam is very small, 
0.808 %, 0.117 %, and 0.646 % for the three LFI bands. 

16.2. Noise and calibration 

Calibration of the TOI (to convert to units of //Kcmb s^^^) 
has improved in several ways since the previous release 
(Planck Collaboration II 2014). Firstly, Planck calibration is 
now based on the dipole signal induced by the annual motion 
of the satellite around the Sun (the orbital dipole). The calibra¬ 
tion thus does not depend on WMAP measurements of the larger 
solar dipole, and it is also absolute, in the sense that it depends 
only on well-measured properties of the solar system and funda¬ 
mental constants. 

Secondly, LFI calibration is now based on full An convolu¬ 
tion of the beam with the dipole (see Sect. 7.1). While the cali¬ 
bration is based on the dipole, the dipole signal is removed from 
the TOI before mapmaking. 

A major source of potential systematic error in calibration 
is Galactic straylight (Galactic emission leaking into the LFI 
horns). We model this effect, and correct the TOI accordingly 
(Sect. 7.4). Straylight (if not corrected) produces evident rings 
centred on the Galactic centre (see Fig. 6). 

As noted, Planck calibration is carried out on a large-scale 
source, namely the orbital dipole, which has a thermal spec¬ 
trum. When assessing the brightness temperature or flux density 
of other astronomical objects with non-thermal spectra, small 
colour corrections are necessary; these are provided in Sect. 7.5. 
For compact sources, the small amount of power missing from 
the main beams, listed above, must be taken into account. As 
an example, the flux density of a compact source with spectral 
index -0.5 extracted from the 30 GHz map requires a 1.00808 


multiplicative correction for missing power and a multiplicative 
colour correction of 0.997. 

The noise properties (white noise levels and knee frequen¬ 
cies) of the LFI receivers are discussed in Sect. 8. The white 
noise was stable over the four-year mission for all receivers. 

16.3. Maps 

LFI produces full-sky maps in Stokes parameters /, Q, and U at 
all three frequencies; the map properties are listed in Sect. 10. 
Calibrated TOI data are destriped using the Madam mapmak¬ 
ing code and maps are constructed using the same package (see 
Sect. 9 for a description and Planck Collaboration VI (2016) for 
full details). In destriping, a mask is employed to limit noise in¬ 
troduced by Galactic emission; the final maps, however, cover 
the entire sky (at Aside = 1024 resolution). Madam also produces 
the noise covariance matrix (NCVM) for each pixel of the maps. 

We also provide maps at lower resolution (Aside = 16; Fig. 9) 
for use in the construction of the low-7’ likelihood (fully de¬ 
scribed in Planck Collaboration XI (2016)). The downgrading 
scheme to smooth the maps from Aside = 1024 to 32 and then 
to 16 is described in Sect. 9.1 and 9.2. Sect. 9.3 describes the 
NCVM for these low-resolution products. 

16.4. Polarization 

The major new feature of this release is the set of polarized 
maps and products. The low-resolution polarization maps at 
70 GHz, in particular, play a crucial role in the construction 
of the Planck low-7 likelihood (Planck Collaboration XI 2016) 
and consequently on Planck values for cosmological parameters 
(Planck Collaboration XIII2016). We therefore devote consider¬ 
able attention to investigating potential systematic errors in these 
maps (detailed in Sect. 11). The largest source of uncertainty 
in LFI polarization measurements is leakage from temperature 
to polarization. This leakage is largely caused by differences in 
the frequency responses or bandpasses between the two arms of 
a given LFI radiometer (“bandpass mismatch”). This mismatch 
can be quantified by a single parameter; for the 70 GHz radiome¬ 
ters, it varies between 0.18 and 1.24 %. The bandpass mismatch 
correction maps are provided in these release at Aside = 256, 
those should be applied to LFI Q and U maps. 

16.5. Validation 

We employ suites of both null tests and simulations to assess the 
quality of LFI maps and other products derived from them (see 
Sect. 12). The null tests exploit the many ways in which the data 
can be divided: survey by survey, year by year, and on the much 
shorter time scale of half-ring differences. The results of some 
of these null tests are shown in Fig. 15; further details appear 
in (Planck Collaboration III 2016) and (Planck Collaboration V 
2016). We call attention to the substantially lower residuals (and 
cleaner maps) resulting mainly from better calibration. The null 
tests do, however, reveal larger than average residual signals in 
polarized maps made from Survey 2 and Survey 4 data (we re¬ 
turn to this issue below). 

Another type of null test is to compare the CMB power spec¬ 
tra derived from different frequencies. This topic is discussed for 
the entire mission in (Planck Collaboration I 2016). Here, we 
point out that Fig. 17, shows good agreement among the three 
LFI bands. In the 7 range 40-300 (which covers the first peak 
of the CMB power spectrum), the three LFI power spectra agree 
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to better than 1 % This agreement extends to measurements of 
compact sources (which involve both a wider I range and val¬ 
ues for the beam solid angles; see Planck Collaboration XXVI 
2016). 

We validated LFI polarization maps by comparing our 
bandpass-mismatch-corrected maps to maps constructed using 
the 72procedure (Sect. 11.4), and found that the Stokes Q 
and U maps were indistinguishable. In particular, the polarized 
structure along the Galactic plane is reproduced, including the 
most significant discrepancies with WMAP maps. 

Simulations based on FFP8 (Planck Collaboration XII2016) 
are also used to validate LFI results. We perform end-to-end sim¬ 
ulations primarily to test the impact of systematic errors and var¬ 
ious steps in our calibration and mapmaking procedures on the 
final results. Sect. 12.6 and Tables 18, 19, and 20 summarize 
the sources of systematic error and their effects on LFI maps. 
The far sidelobes of the LFI beams are the dominant source of 
uncertainty in the 30 GHz maps. At 44 and 70 GHz, other instru¬ 
mental effects dominate, particularly 1 Hz electronic spikes and 
ADC non-linearity, respectively. The overall systematic effects 
uncertainty was estimated to be 0.88, 1.97, and 1.87 /iKcMBin 
the 1 component; 1.11, 1.14, and 2.25 yuKcMBin the Q compo¬ 
nent; 0.95, 1.20, and 2.22 in the U component at 30, 44 and 
70 GHz, respectively. 

As mentioned above, null tests show that the polarized data 
from Surveys 2 and 4 contain residual signals (possibly due to 
contamination from Galactic emission). As a consequence, we 
choose to be conservative and omit these two surveys (approx¬ 
imately 1/4 of the data) from the low-7’ likelihood. The studies 
supporting this decision are described in Sect. 13. The optical 
depth, T, is the cosmological parameter most affected; including 
or omitting Surveys 2 and 4 changes r by about 0.5 cr. 
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